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FOREWORD 


The  Takeoff  and  Landing  Digital  Computer  Program  was  prepared  by  the 
Convair  Aerospace  Division  of  General  Dynamics  Corporation  under  USAF 
Contract  F33615-71-C-1754,  Project  643A,  "STOL  Tactical  Aircraft 
Investigation. "  This  contract  was  sponsored  by  the  Prototype  Division  of 
the  Air  Force  Flight  Dynamics  Laboratory.  The  USAF  Project  Engineer 
was  G.  Oates  (PT)  and  the  Convair  Aerospace  Program  Manager  was 
J.  Hebert.  C.  A.  Whitney  was  the  principal  contributor. 

The  research  reported  was  conducted  during  the  period  from  7  June  1971 
through  31  January  1973.  This  report  was  submitted  by  the  author  on 
31  January  1973  under  contractor  report  number  GDCA-DHG73-001. 

This  report  has  been  reviewed  and  is  approved. 
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ABSTRACT 


The  MILSTOL  (MILitary  STOL)  takeoff  and  landing  digital  computer  program  was 
developed  under  USAF  Contract  F33615-71-C-1754,  "STOL  Tactical  Aircraft  Investi¬ 
gation,  "  to  compute  takeoff  and  landing  characteristics  of  powered-lift  STOL  aircraft. 
S  calculates  a  point  mass  takeoff  and/or  landing  for  a  trimmed  configuration  with 
either  externally  blown  jet  flaps,  internally  blown  jet  flaps,  or  mechanical  flaps  with 
vectored  thrust  within  the  constraints  set  forth  in  Reference  1.  Contained  in  this 
report  are: 

1.  Discussion  of  assumptions  and  methods  used  in  the  trajectory  calculations. 


2,  Definition  of  common  list  variables. 


3.  Definition  of  the  input  variables  and  sample  input  data  for  the  externally 
blown  jet  flap  configuration. 

4.  Sample  output  for  the  externally  blown  flap  configuration. 


5.  Program  listings  and  flow  charts. 

The  program  is  wrftLa  in  Fortran  IV  for  use  on  CDC  6000  series  digital  computers 
,u.,  requires  37  K  central  memory  for  loading  and  execution.  Hie  program  is  com- 
patable  with  both  the  CDC  RUN  and  FTN  compiler  systems. 
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SECTION  1 
INTRODUCTION 


The  military  STOI.  takeoff  and  landing  digital  computer  program  (MILSTOL  -  Convair 
Aerospace  Division,  San  Diego  Operation,  scientific  computer  program  P6592)  was 
developed  to  calculate  the  takeoff  and  landing  performance  of  powered-llft  STOL  air¬ 
craft.  The  performance  calculations  are  made  using  exact  two-degree-of-freedom 
equations  of  motion  for  a  point  mass  aircraft  (i.e. ,  no  pitch  dynamics).  The  program 
was  developed  from  the  takeoff  portion  of  the  Aircraft  Performance  Analysis  System, 
Reference  2,  and  uses  data  handling,  equations  of  motion,  and  general  use  subroutines 
from  that  program. 


R  The  takeoff  portion  of  the  MILS’!  OL  program  performs  a  cor  stant-weight  "balanced" 

takeoff  from  zero  forward  speed  to  liftoff  and  to  stop,  within  the  constraints  shown  in 
R  Figure  1-1,  for  a  matrix  of  gross  weights  and  runway  altitudes.  Velocity  cues  for 

§  engine  failure,  rotation,  and  liftoff  are  factors  times  the  minimum  control  speed  and 

I  the  stall  speed  with  power  on  and  the  critical  engine  failed. 


The  landing  phase  of  this  program  performs  a  "no-flare"  style  landing  approach, 
touchdown,  and  deceleration  to  stop  for  the  input  matrix  of  gross  weights  and  runway 
altitudes.  Velocity  cues  for  approach  and  touchdown  speeds  are  functions  of  the 
minimum  control  speed  and  power-on  stall  speed  with  the  critical  engine  failed.  The 
landing  trajectory  is  calculated  within  the  constraints  shown  in  Figure  1-2.  Aero- 
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Figure  1-1.  Balanced  Field  Takeoff  Ground  Rules  for  STAI* 
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dynamic  and  propulsion  data  is  input  in  tabular  form  and  is  bandied  by  individual 
modularised  subroutines.  The  atmospheric  properties  subroutine  is  compatible  with 
die  1962  U<  S.  Standard  Atmosphere  and  the  MIL-STD--210A  temperature  conditions. 
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Figure  1-2.  Landing  Ground  Bales  lor  STAL 
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SECTION  2 

PROGRAM  DISCUSSION 


The  subroutines  in  the  MILSTOL  program  are  classified  into  five  functional  cate¬ 
gories. 

].,  Executive  Program 

2.  Maneuver  Driving  Subroutines 

3.  Physical  Data  Subroutines 

4.  General  Use  Subroutines 

5.  Data  Handling  Subroutines 


Program  flow  and  structure  are  shown  in  Figure  2-1;  each  program  subroutine  is 
discussed  by  functional  category  in  the  following  sections. 


Figure  2-1.  MILSTOL  Functional  Structure 
2-1 


2.  i  EXECUTIVE  PROGRAM  (MILSTOL) 


The  Executive  Program  controls  the  reading  of  inputs,  the  initializing  of  weights, 
altitudes,  and  temperature,  and  execution  of  the  appropriate  trajectory  subroutine 
for  the  given  weight  and  altitude  matrix.  In  addition,  a  Dump  Overlay  program  is 
included  in  the  MILSTOL  procedure.  This  Dump  Overlay  program  is  executed  only 
in  the  event  of  a  fatal  error.  At  the  time  of  the  abnormal  termination,  the  Dump 
Overlay  program  is  loaded  from  a  local  file  and  prints  a  listing  of  all  common  list 
variables  in  Namelist  format. 

2.2  MANEUVER  DRIVING  SUBROUTINES 

Two  maneuver  subroutines,  TAKEOFF  and  LANDING,  are  incorporated  in  the 
MILSTOL  program. 

2. 2. 1  TAKEOFF  TT  "TORY  SUBROUTINE  (TAKEOFF)  -  The  TAKEOFF  sub¬ 
routine  is  the  driver  for  all  portions  of  the  takeoff  maneuver.  Ground  rules  and 
constraints  for  this  maneuver  are  shown  in  Figure  1-1.  This  subroutine  reads  tra¬ 
jectory-related  variables  (e.g. ,  rolling  and  braking  coefficients  of  friction,  time  de¬ 
lays,  minimum  control  speeds,  etc. )  executes  the  takeoff  calculation  procedure,  ad¬ 
justs  the  engine  failure  speed  to  balance  the  continued  and  aborted  takeoffs,  and  causes 
the  takeoff  time  history  to  be  output. 

Because  of  the  balanced-type  takeoff  specified  by  the  ground  rules  and  to  ensure 
efficient  program  operation,  the  takeoff  is  not  calculated  as  a  continuous  function  of 
time.  Instead,  the  program  is  divided  into  segments.  The  sequence  of  calculations 
is: 

1.  Stall  speed  with  power  on  and  the  critical  engine  failed  is  calculated.  Liftoff 
speed  is  set  using  the  air  minimum  control  speed  and  stall  speed.  The  initial 
value  of  engine  failure  speed  is  set  equal  to  ground  minimum  control  speed. 

2.  Angle  of  attack  for  liftoff  is  calculated  with  the  critical  engine  failed.  If  the 
ground  contact  angle  (the  angle  for  the  tail  striking  the  ground  during  rotation) 
is  exceeded,  liftoff  speed  is  increased  by  one  percent  of  the  air  minimum 
control  speed  or  power-on  stall  speed.  At  this  point,  the  pertinent  aircraft 
conditions  are  output.  When  the  conditions  for  liftoff  are  established,  maxi¬ 
mum  rate  of  climb  at  the  liftoff  speed  is  calculated  and  output, 

3.  The  first  segment  of  the  takeoff  is  calculated  with  the  critical  engine  foiled 
by  integrating  time  and  tangent' al  acceleration,  along  with  the  input  rotation 
rate  in  a  negative  sense  from  liftoff  speed  to  rotation  speed  to  obtain  velocity, 
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distance,  and  aircraft  attitude.  The  integration  is  terminated  when  the  air¬ 
craft  pitch  attitude  is  zero,  and  a  check  is  made  to  ensure  that  the  rotation 
velocity  is  greater  than  the  selected  engine  failure  speed  and  the  ground  mini¬ 
mum  control  speed.  3f  this  criterion  is  not  satisfied,  the  liftoff  speed  is 
increased  using  an  empirical  relationship.  The  program  then  returns  to 
Step  2  and  continues  until  the  rotation  velocity  criterion  is  satisfied. 

4.  If  rotation  velocity  is  greater  than  engine  failure  speed,  time  and  tangential 
acceleration  are  integrated,  in  a  negative  sense  from  the  rotation 
velocity  to  the  engine  failure  speed,  with  the  critical  engine  failed  to  obtain 
velocity  and  distance.  This  distance  is  the  second  takeoff  segment  and,  when 
added  to  the  first  segment  from  Step  3,  is  the  "continued  takeoff"  distanoe 
used  in  the  field-balancing  relationship. 

5.  The  distance  for  the  "aborted  takeoff"  used  in  balancing  the  field  length  is 
calculated  by  integrating,  in  a  positive  sense,  the  time  and  tangential  ac¬ 
celeration  valuables  from  engine  failure  speed  to  stop  for  velocity  and  dis¬ 
tance  with  the  critical  engine  failed.  During  the  integration,  engines  are 
set  to  idle  after  the  reaction  time  and  the  deceleration  devices  (brakes,  lift 
dumpers,  and  reverse  thrust)  are  deployed  at  the  end  of  an  actuation  time 
interval. 

6.  The  "continued  takeoff"  distance,  Steps  3  and  4,  and  the  "aborted  takeoff" 
distanoe  are  then  used  in  a  linear  convergence  procedure  to  adjust  the  engine 
failure  speed  so  that  these  distances  are  equal.  After  the  new  engine -failure 
speed  is  selected,  the  program  returns  to  Step  2  and  calculates  new  "continued" 
and  "aborted"  takeoff  distances.  If  the  "aborted  takeoff"  distance  is  greater 
than  or  equal  to  the  "continued  takeoff"  distanoe  and  the  engine  failure  speed 

is  equal  to  the  ground  minimum  control  speed,  the  takeoff  is  by  definition 
balanced  and  the  program  continues. 

7.  After  the  preceding  stops  have  balanced  the  takeoff  distance,  time  and 
tangential  acceleration  are  integrated  from  start  to  the  engine 

failure  speed  (with  all  engines  operating)  for  velocity  and  distanoe.  Sum¬ 
mations  of  the  distances  from  Steps  3  and  4  with  this  distanoe  and  the 
distanoe  from  Step  5  comprises  the  balanced  takeoff  distanoe  for  this  con¬ 
figuration. 

2.2.2  LANDING  TRAJECTORY  SUBROUTINE  (LANDING)  -  The  LANDING  sub¬ 
routine  is  the  driver  for  all  portions  of  the  landing  maneuver.  Ground  rules  and 
constraints  for  the  landing  trajectory  are  shown  in  Figure  1-2.  This  subroutine  reads 
trajectory-related  inputs,  executes  the  landing  calculation  procedure,  and  outputs  a 
landing  trajectory  time  history  for  the  input  configuration. 
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TLe  sequence  of  calculations  is: 

1.  Stall  speed  with  power  on  and  the  oritioal  engine  failed  is  calculated.  The 
approach  speed  that  conforms  to  the  ground  rules  of  Figure  1-2  is  set  using 
the  air  minimum  control  speed  ami  the  power-on  stall  speed. 

2.  Angle  of  attack  at  touchdown  is  calculated  with  all  engines  operating.  If  the 
pitch  attitude  exceeds  the  ground  contact  angle,  the  approach  speed  is  in¬ 
creased  by  one  percent  of  the  power-on  stall  speed.  If  the  rate  of  sink  at 
touchdown  exceeds  the  input  maximum,  the  glideslope  angle  is  reduced  so 
that  the  rate  of  sink  limit  is  met.  If  the  configuration  attitude  is  such  that  the 
nosewheel  hits  first,  the  program  prints  an  error  message  and  returns 
control  to  tb*  executive  routine.  When  all  touchdown  criteria  are  satisfied, 
the  program  outputs  the  conditions  at  touchdown  and  calculates  the  maximum 
rate  of  climb  available  at  touchdown  speed  with  the  critical  engine  failed. 

U.  Tb~  program  then  calculates  angle  of  attack  at  the  obstacle  (with  all  engines 
operating).  This  calculation  is  performed  at  the  approach  speed  calculated 
in  Step  2.  Aircraft  conditions  at  the  obstacle  are  output  along  with  the  maxi¬ 
mum  rate  of  climb  available  with  the  critical  engine  failed.  This  calculation 
is  made  at  the  obstacle  to  account  for  ground  effects. 

4.  Landing  air  distance  is  calculated  by  performing  a  one-step  integration  using 
the  velocities  from  Steps  2  and  3. 

5.  After  touchdown,  the  program  performs  a  step-wise  integration  from  touch¬ 
down  to  stop  to  calculate  ground  distance.  During  this  integration,  engines 

are  set  to  idle  power,  aircraft  attitude  is  rotated  down  to  zero,  and  deceleration 
devices  (brakes,  lift  dumpers,  and  reverse  thrust)  are  deployed,  after  allow¬ 
ing  for  actuation  time  delay, 

6.  Total  landing  distance  Is  the  summation  of  distances  from  Steps  4  and  5, 

2. 3  PHYSICAL  DATA  SUBROUTINES 

ATMOS,  AEROI,  and  THSTI  are  the  three  physical  data  subroutines  in  the  MILSTOL 
program.  These  subroutines  read  inputs  and  store,  retrieve,  and  calculate  the 
atmospheric,  aerodynamic,  and  propulsion  characteristics  required  to  solve  the 
equations  of  motion. 

2.3.1  ATMOSPHERIC  PROPERTIES  SUBROUTINE  (ATMOS)  -  The  ATMOS  sub¬ 
routine  supplies  the  program  with  ambient  temperature,  pressure,  density  ratio, 
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and  speed  of  sound  as  a  function  of  altitude  and  type  of  day  or  an  input  temperature. 

The  temperature/type  of  day  options  available  are: 

1.  U.S.  Standard  Atmosphere,  1962. 

2.  MIL-STD-210A  Tropic  Day. 

3.  MIL-STD-210A  Polar  Day. 

4.  MIL-3TD-210A  Hot  Day- 

5.  MIL-STD-210A  Cold  Day- 

6.  An  input  temperature  in  °F. 

These  options  all  use  the  standard  day  pressure  altitude  relationship  in  the  cal¬ 
culation  procedure.  Options  1  through  5  use  the  appropriate  temperatures  from 
References  3  and  4. 

2. 3. 2  AERODYNAMIC  DATA  SUBROUTINE  (AEROl)  -  This  subroutine  was  develop¬ 
ed  to  store  and  retrieve  trimmed  aerodynamic  data  for  configurations  with  externally 
blown  jet  flaps,  internally  blown  jet  flaps,  and  mechanical  flaps  with  vectored  thrust. 
Conventional  configurations  without  thrust  augmented  lift  can  be  used  by  either  modify¬ 
ing  the  table  lookup  procedure  or  by  entering  the  power-off  data  for  four  dummy  thrust 
coefficients  and  using  the  mechanical  flap  plus  vectored  thrust  option. 

AEROl  has  two  entries:  AEROl ,  which  retrieves  maximum  lift  characteristics  as  a 
function  of  flap  deflection  and  momentum  coefficient,  and  AER02,  which  retrieves 
lift  and  drag  data  as  a  function  of  flap  deflection,  angle  of  attack,  and  momentum 
coefficient.  The  inputs  to  this  subroutine  are  geometrio  data,  configuration  type, 
and  aerodynamio  data  tables. 

Because  of  the  differences  in  methods  for  estimating  aerodynamic  data  for  each  of 
the  three  configurations,  there  is  a  unique  method  of  storing  and  retrieving  the  data 
for  eaoh  of  the  three.  The  externally  blown  jet  flap  data  indudes  all  direct  and  in¬ 
direct  thrust  effects.  The  internally  blown  Jet  flap  data  includes  all  thrust  effects 
due  to  trailing  edge  slot  blowing  but  none  of  the  thrust  effeots  duo  to  the  cruise 
engines.  The  mechanical  flap  plus  vectored  thrust  sequenoe  assumes  that  the  aero¬ 
dynamic  data  includes  all  indirect  thrust  effects  including  any  supercirculation 
effects,  but  none  of  the  direot  thrust  vector  effects. 

All  additional  items  that  degrade  the  aerodynamio  data  (e.g. ,  lift  dumpers, 


2-5 


engine  out  corrections,  etc.)  are  cued  from  the  trajectory  subroutines  and  are 
included  in  the  final  lift  and  drag  values  before  returning  to  the  calling  subroutine. 

Retrieval  and  calculation  of  the  aerodynamic  characteristics  with  the  critioal  engine 
failed  are  handled  in  the  same  manner  for  all  configurations.  In  this  calculation 
scheme,  it  is  assumed  that  the  input  aerodynamic  data  tables  are  valid  for  either 
the  all-engines -operating  or  the  one -engine -failed  condition  if  a  later  correction 
is  made  to  compensate  for  engine-out  moments.  The  correction  for  trimming  of 
engine-out  moments  are  input  increments  to  lift  and  drag.  Configuration  design 
and  aerodynamic  conditions  that  allow  this  assumption  are: 

1.  All  internally  blown  jet  flap  configurations  are  assumed  to  have  cross-ducting 
so  that  the  spanwise  distribution  of  blowing  is  always  symmetric. 

2.  Engine-out  degradation  in  the  mechanical  flap  plus  vectored  thrust  con¬ 
figurations  is  due  to  the  loss  of  the  thrust  vector  component,  not  loss  of 
aerodynamic  lift. 

3.  For  configurations  with  externally  blown  jet  flaps,  the  aerodynamic  charac¬ 
teristics  are  a  function  of  total  thrust  coefficient.  Analysis  of  the  data  in 
Reference  5  shows  that  aerodynamic  lift  and  drag  characteristics  are  the 
same  for  a  given  total  thrust  coefficient,  whether  or  not  an  engine  is  failed, 
if  the  moments  due  to  the  engine  failure  are  not  trimmed.  The  lift  and  drag 
increments  due  to  trimming  the  engine-out  rolling  and  yawing  moments  are 
incorporated  as  stated  above, 

2.3.3  PROPULSION  CHARACTERISTICS  SUBROUTINE  (THST1)  -  The  THST1  sub¬ 
routine  is  used  for  storing  and  retrieving  single -altitude  propulsion  characteristics 
for  all  three  powered  lift  configurations.  Inputs  to  this  routine  are  in  the  form  of 
single -altitude  maximum  gross  thrust,  ram  drag  at  maximum  thrust,  idle  gross 
thrust,  maximum  reverse  thrust,  windmilling  drag  of  a  single  engine,  and  maximum 
gross  thrust  at  the  slot  exit  for  internally  blown  jet  flap  configurations  as  a  function 
of  velocity.  Procedures  have  been  incorporated  so  that  reverse  thrust  is  always 
symmetric  in  the  engine-out  case. 

2.4  GENERAL  USE  SUBROUTINES 

There  are  eight  general  use  subroutines  in  the  MILSTOL  program.  Those  are 
general  purpose  subroutines  extracted  from  Reference  2  and  have  applications  outside 
of  the  MILSTOL  program.  These  subroutines  perform  integrations,  print  page 
headings,  handle  equations  of  motion,  and  provide  the  logic  for  convergence  pro¬ 
cedures,  Each  of  these  subroutines  is  described  in  alphabetical  order  in  the  follow¬ 
ing  paragraphs. 
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2.4. 1  AEQFM  —  To  maintain  consistency  between  different  programs  and  cal¬ 
culations,  all  equation-of-motion  calculations  are  performed  using  the  AEQFM 
subroutine.  Figure  2-2  shows  the  axis  systems,  foroes,  and  angles  used  in  the  two- 
degree-of -freedom  calculations.  Equations  1  and  2  are  balanced  for  the  appropriate 
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Figure  2-2,  Force  and  Angle  System  Used  In  the  MILSTOL  Program. 

flight  condition.  When  acceleration  or  deceleration  on  the  ground  is  required,  an 
additional  term  is  added  to  Equation  1  to  account  for  the  ground  friction  force.  All 
accelerations  in  Equations  1  and  2  arc  converged  simultaneously  using  the  TRIM  and 
SMLT2  subroutines. 

2,4, 2  GILL  —  This  is  an  Integrating  subroutine  that  uses  the  method  developed  by 
S.  Gill  (Reference  6),  to  provide  fourth -order  accuracy  while  requiring  a  minimum 
number  of  storage  registers.  The  subroutine  requires  four  passes  to  accomplish 
the  integrating  step; 

Pass  One.  Take  derivatives  at  the  start  of  the  interval  and  predict  conditions 
at  the  middle  of  the  interval. 

Pass  Two,  Take  derivatives  based  ou  predicted  conditions  at  the  mid-interval 
and  combine  with  derivatives  from  first  pass  to  predict  conditions  at  the  mid¬ 
interval. 
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Pass  Three.  Take  derivatives  based  on  latest  estimate  of  mid-interval 
oonditions  and  combine  with  derivatives  from  first  two  passes  to  predict  end-of- 
interval  conditions. 

Pass  Four.  Take  derivatives  based  on  ent*  'f-interval  conditions  and  combine 
with  derivatives  from  other  passes  to  calculate  conditions  at  the  end  of  the 
interval. 

This  process  is  repeated  for  each  integration  step.  The  calling  subroutine,  GINTG, 
keeps  track  of  the  number  of  passes  and  checks  for  terminations  after  four  passes. 

2. 4. 3  GINTG  —  This  subroutine  is  the  driver  for  trajectory  integrations.  It  calls 
the  equations -of-motion  subroutine  (AEQFM)  for  accelerations,  then  the  GILL  sub¬ 
routine  to  integrate  for  velocities,  distances,  and  aircraft  attitude. 

2. 4.4  HEAD  SUBROUTINE  —  This  subroutine  is  used  for  printing  columnar  headings 
in  the  time  history  printout. 

2.4.5  NWRP2  —  This  is  a  Newton- Wrapson  iteration  subroutine  which  determines  the 
value  of  x  that  will  return  y  equal  to  zero  based  on  a  linear  prediction  using  two 
previously  calculated  points. 

2. 4. 6  SKIP  —  This  is  an  output  formatting  subroutine  that  starts  a  new  page,  prints 
a  standard  page  heading,  and  restarts  the  line  count. 

2.4.7  SMLT2  —  This  subroutine  performs  a  simultaneous  equation  solution,  using 
derivatives  from  subroutine  TRIM,  to  obtain  increments  to  the  independent  variables 
that  will  result  in  the  desired  accelerations  for  subroutine  AEQFM. 

2. 4. 8  TRIM  —  This  subroutine  controls  tho  systematic  perturbation  of  independent 
variables  and  stores  the  variation  of  each  acceleration  with  respect  to  each  variable 
as  a  derivative. 

2. 5  DATA  HANDLING  SUBROUTINES 

The  MILSTOL  program  usee  five  subroutines  and  two  functions  whose  solo  purpose 
is  data  handling.  These  routines  fit  curves  to  data,  evaluate  curve  fits,  and  per¬ 
form  tabic  look-upe.  These  routines  and  functions  are  discussed  in  alphabetical 
order  in  the  following  paragraphs. 

2. 5. 1  ALIN  —  This  function  makes  a  linear  fit  between  two  ( x,y)  points. 


2. 3. 2  FIND  —  This  subroutine  performs  a  one-dimensional  linear  interpolation 
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within  a  data  table.  For  arguments  greater  than  or  less  than  the  table,  a  linear 
extrapolation  is  performed. 

2. 5. 3  INTP2  -  Subroutine  INTP2  fits  a  third-order  polynomial  through  four  (x»y) 
points,  returns  the  coefficients,  and  returns  a  y-answer  for  an  x-argument. 

2. 5.4  KABD  -  The  subroutine  evaluates  a  hyperbolic  fit  by  returning  a  y-answer 
for  an  x-argument  of  the  equation  y  =  K/(x-A)  +  B  +  D  x,  where  coefficients,  K,  A, 

B,  anl  D  are  provided  in  the  calling  list. 

2. 5. 5  LAGRA  —  Subroutine  LAGRA  returns  a  y-answer  for  an  x-argument  using 
a  Lagrangian  interpolation  on  four  (x,y)  points. 

2. 5. 6  LOOK  —  Originally  written  to  handle  three-dimensional  tabulated  thrust  and 
fuel  flow  data,  this  subroutine  has  been  developed  into  a  more  general  form  that  can 
handle  any  three-dimensional  tabular  data.  The  LOOK  subroutine  performs  the  table 
lookup  using  a  non-linear  technique  the  basis  of  which  is  the  LAGRA  subroutine.  S 
has  four  options  for  locating  data  and  it  can  also  return  derivatives  with  respect  to 
three  independent  variables  using  a  four-point  interpolation  of  each  independent  vari¬ 
able. 

2. 5.7  SL  —  Ihe  SL  function  calculates  the  linear  slope  between  two  (x,y)  points. 
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COMMON  LIST  VARIABLES 

Tuia  section  contains  three  tables  that  provide  the  user  with  a  key  to  the  definition 
and  usage  of  the  labeled  common  blocks  incorporated  in  the  MILSTOL  program. 
Table  3-1  describes  each  labeled  common  block.  The  blocks  have  been  constructed 
by  function  to  aid  the  user  in  future  modifications  or  upgrading. 

Table  3-1.  Description  of  Labeled  Common  Blocks. 

LIST  1  Contains  variables  used  for  input  and  output  units,  carriage  control, 
and  page  headings. 

LIST  2  Contains  variables  describing  forces,  velocities,  altitudes,  and  weights, 

LIST  3  Contains  coefficients ,  angles,  and  aerodynamically  significant  geometry. 

LIST  4  Contains  variables  used  for  transmitting  propulsion  characteristics, 

LIST  5  Contains  variables  used  for  atmospheric  properties. 

LIST  6  Contains  physical  constants  and  conversion  factors. 

LIST  9  Used  for  transmitting  data  from  the  LOCK  subroutine. 

LIST  15  Contains  integration  variables  and  controls. 

LGEOM  Contains  aircraft  geometry  and  angles, 

CONTROL  Control  flags  used  duri  ^  takeoffs  and  landings, 

LEST  99  Error  index  Sag. 
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Table  3-2.  Subroutines  In  Which  Labeled  Common  Blocks  Are  Used. 


Common  List  Name 


g 

L  IST  i 

LIST  2 

LISTS 

LIST  4 

LISTS 

to 

B 

a 

LIST  8 

LIST  9 

LIST  15 

LGEOM 

o 

o 

LIST  99 

MILSTOL 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

(Main  Program) 

MILSTOL 
(Dump  Program) 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

AEQFM 

X 

X 

X 

X 

X 

X 

X 

X 

X 

AEROl 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

ATMOS 

X 

X 

X 

FIND 

X 

GINTG 

X 

X 

X 

X 

X 

X 

X 

X 

X 

HEAD 

X 

KABP 

X 

LANDING 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

LOOK 

X 

SKIP 

X 

X 

X 

X 

TAKEOFF 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

TH5TI 

X 

X 

X 

X 

X 

X 

X 

X 
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Table  2-3.  liatinitlsa  of  Common  Block  Variables. 
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V*OUPCC 

CO^wOV 

oepinition  \  */a 

NAV# 

NAME 

ALrap 

LIST* 

WING  ANGL -  OP  ATTACH  (OADIANJt 

AU*»X© 

USTJ 

HA* I Mu  .  ACCCpASlE  ANGLE  OP  ATTACH  (MCI  ANSI 

At*>HC 

CIST3 

»|KO  ANGLE  LP  ATTACK  poo  STALL  IPAOJAKS;  . 

AC  TWO 

LIST3 

ANGLE  OP  ATTACH  UP  The  tM'ilne.  -Lie  MCA43IAMOI 

iw 

lists 

Af.CeLSBATICW  itCKKAL  TO  THE  VELOCITt  VCCTQA  IPT/SCC  401 

Awccr 

CGfC* 

H'J-r  USSO 

AK? 

CI5TO 

the  value  or  »e»  betusnfd  pgom  sutniouiiNe  coon 

A*<2 

LIST" 

inf  value  ~r  isfi  oetubnco  reep  sumoctike  coo* 

*0 

LISTS 

wot  used 

AT 

cute 

Arcci.'-OKTiec  pap*,.lEc  to  the  vEcocftv  vecroa  ipt/sec  soi 

IX 

•  isre 

ACCTuTPATtCN  r-AOALLEC  TC  THE  EABSN  IP T/SCC  SOI 

1? 

L !  STO 

ACCELERATION  W»L  TO  IMC  EABTH  tPTESCC  SOI 

A* 

CISTS 

WING  SPAN  (PTI 

CAOf> 

USTI 

IINTEGEAI  UNIT  POP  nCAOING  CAPO  INPUT 

CPAO 

CGIO'I 

NOT  CISCO 

€C* 

usri 

OOAO  COEFFICIENT 

C**CC 

LISTS 

NOT  USED 

cu 

Cl  ST 3 

CIPT  COfPf  ICIENT 

CISTS 

NOT  WEED 

ccr* 

CISTS 

COEFFICIENT  OP  EVICTION 

Ct  AO 

CISTS 

NOT  USED 

Ct^AH 

CISTS 

HA* 1*0*  CIPT  COEPP ICIENT 

cc«*c 

cc^c4- 

nct  useo 

CC5TW 

1.6TO'* 

NOT  USE 9 

c* 

LISTS 

not  urea 

cr 

CUTS 

NOT  USED 

c  AWT 

LlStiS 

height  op  ming  AfOve  the  wound  im 

OFT* 

LISTS 

NOT  us  Q 

p*rv 

CISTS 

NOT  USED 

t^r  - 

LISTS 

nc*  useo 

LISTS 

not  u*eo 

Wf*<‘ 

CISTS 

not  used 

LISTS 

nST  USCO 

O^tTO 

L>l  " 

not  UfEO 

O*\o 

LI ST  Is 

r>.jo  opplsction  idegi 

C-C** 

LIST? 

».-  — <  -K4V1C  rcecr  IN  THt  9AIA4  v  >  PCCT  ICn  IcOSI 

6**w' 

L  ISTS 

not  u*fr 

LI»TS 

ti«T  USP* 

D*tw* 

LOCO" 

nct  ussa 

CT 

LIST? 

NOT  L-SSO 

LETCH 

not  uSSo 

CT**>W 

c«c* 

NOT  LrStO 

ft***# 

NO*  uf»*> 

L 1 3 1 1  * 

*  1 VT  INTfiOtTlNG  INTCnvSL  l«^tl 

Q.**« 

LISTS 

EAnv*3Ts  ?*C*r? t  *e  a^to-'S 

£»• 

LIST* 

ar»  ?p*5  pee  JLV  T*etiLET-  Ac»*n  ?“t  vEt??'**  vCETO*  iljpsi 

8*V*/ 

C«0* 

NET  u?.?o 

f*.C*£ 

CUT* 

*«p  NV'flfP  OP  fNElNTT  L'PvfcSTINC 

xr« 

CIST* 

not  v<5#a 

i*r 

CISTf» 

fstTft  »(K?f  IN  T»f  CIPT  lieftllaN  fCHSl 

*Tr* 

LISTS 

«)Nv*e*t  f«f/Sft*A5  *0  MN5TJ 

*TfW 

CISTS 

SCNVE<i«t  PTtt  TO  NtuTICtL  PIC** 

L  1ST? 

N®S  »1T0 

L  l&*  1* 

ANrtLt  <SP  ATTACK  PCS*  iESNPrf'  CONTACT  ILt$||L|l 

£**«* 

CIST.' 

♦light  Path  pel!  ‘»AO|An-»i 

ft** 

C«S»t 

L;a*  Pi<»0»  v-r-TAL  TO  T«e  ViCCCIf*  VtCTOP  tOl 

ft* 

lists 

L«*J  TACTCS  PAi.'ALVSC  *«  *“i  VT.L»i|*»  VfiTO*  let 

ft? 

CIST* 

4CrTC{t»t«e*  5A.P  TO  V^ASIT*  |PT/J«e  SOI 

V^l 

CI1*I 

VAVIAts.*  NASH  AVAIL A*Cf  PC*  ewTOur  Nisivst 

CU»f 

*L**TuC;t  fp»l 

UJT? 

NCT  W1C0 

in»»» 

LI*T|* 

|N?T*  ?.»»  »**<.  'IM5CS  CP  attach  iS  e«AO«ff 

LIST* 

f«ffl  *9*e 

1  AT 

CSV*  IS 

NT".A*t it  ♦t<*v**A*ioN!t»s:T»  av6i«e  tmaw 

i,n« 

C8**l 

PAT*  l•«^T»4A*'/Nf  **| 

rw 

CtCTfr* 

rsrsj*  t*uNT* 

l*V*#- 

COST* 

not 

u  r*** 

CISTS* 

lM»  i^fA  TV  hA«CVI|Ut  c ICC 

ClST.T 

«4l  Irtto 

t** 

MATE 

tlittTS  JHSIIS  SnieifcjAittNs  TO  PJ A?  OP.  t»S4«TE 

list* 

OPTIC".  AfvCCSOi*  »c«  sv*t-i»urTINE  4COP* 

♦xrfrtft 

C!tT3 

Is!  CP  PUL®  TNf.ie**t»  P'*9  tLfiilNX.'iNE  L to* 

CiS»IS 

SHfCTA  VUU'-LTI  16  lvl*Si*«t  IN  MACVtlNl  Slvtt 

»** 

CIST, 

W*6«luf  call!  e^r  *l!»  rtJiC r  1  - 

ili«l 

NOT  capo 

!«*v 

tcwcwrc 

•IVIAT*  »«ty»T  IMSICATOP 

$*«» 

Ll»«  av"Af4  l«Ot*Al6l 

}**• 

CISTf 

Sl««tvs  CPT»o»  is  aeopn 

cafe* 

nc*  veto 

e«*Ts*c 

iNfiir*  fjo  T»*t  op  *c«f<t6  tip*  tnn* 
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Table  3-3.  Definition  of  Common  Block  Variables.  (Contd) 


w: 


1 


t 


JPTV 

LIST* 

MOT  US tO 

KTCf 

LISTS 

ISE*L)  CCNVEBTS  KNOTS  TO  TC  IT /at  COM3 

LIFT 

LIST2 

IRC4LI  1CBOOVN4KIC  TCSCt  |N  ThE.  Lit?  OlSeCTIQN  Util 

LI«1T 

LISTl 

L  INC  LIMIT  RCP  CuTOuT  P»0e 

Lt*e 

LlSTl 

ThC  CuCBCnT  Of  LINES  vSITtEw  ON  Th*T  paoc 

**CH 

LIST* 

IDEAL!  Tout  KNCM  MU**SCO 

NO 

LISTS 

OPTION  StLCCTOO  ton  SU9B0UTINC  LOOK 

LISTS 

NOT  USEO 

NTM“. 

contsce. 

index  too  BtvTosr  thlust  *oce 

LISTS 

ib£*li  convcbts  nnuticnl  wiles  to  rttt 

Mj 

LISTS 

ioc»Li  kinematic  viscos i Tv  eso.rT/seci 

LIST! 

(INTECEBl  CUOOENT  0»ct  NUKCII.B 

lITTS 

*t«os«n*£oic  *«sieNT  pcessuat  iL9S/60.tTi 

L«eo'1 

**ot  useo 

POINT 

LtSTl 

MMfMOt  VJN|T  tea  .SITING  POlNTtO  OUTPUT 

LISTS 

se*l  level  *T«*oiC)*eoie  pcessust  tLt/so.rti 

0 

LISTS 

*iocp*»'.  pitch  bite  ibaouns/scci 

QOCT 

LCtO' 

not  USED 

0? 

LISTS 

dynamic  popssuot  ti»£s  serestNce  aoen  itesi 

0*C 

LISTS 

NOT  UMO 

flrr 

LIST? 

0»T£  or  CLIP*  ItT/SfC* 

CH407 

lists 

se*  level  5T**c*ao  o*v  nib  cc.Nt.tTv  tSfeU&iscuaic  rti 

0*022 

LISTS 

On02  01 VIDEO  ev  t 

0*7 

lcco* 

NOT  used 

cth 

LGEO-* 

wot  USEO 

orcr* 

lists 

CCNVEBTS  QUO  I  Ann  TO  DECREES 

Ova 

lceo* 

MOT  USED 

s 

lists 

•  INC  LEttaeNCE  Not*  ISO.PTI 

SCO 

lcso* 

NOT  US»0 

$►*«*» 

LOTO'* 

mot  useo 

516 

list; 

DENSITY  BATIO 

S!w*L 

LOSS'* 

NOT  U*SD 

5 

LOT  O'* 

MOT  u*LO 

t  IS*S 

s«*»EO  -or  souma  ipt/seci 

Tfvr 

LISTS 

*|E  Te*u>«a*TuCt  iota  r*xCi:F**EIT) 

LISTS 

.to  »erBEu*TucE  toco  attxtnCt 

tGOO<^ 

LIST* 

on  ss  tmsust  roe  *ll  engines  nlCwo  yhe  tuous*  *»i» 

Tu-rtO 

lots- 

»|OCOIUT  STTJTUOC  IMJIIKSI 

tM|» 

L  ISTS 

tmoust  iweiafwTf  i3*oi«n$i 

LIST* 

CBOSt  thOuSt  »t  ynt  just  r»iT  rc*»  I*r  cowries  iu»S> 

LIST* 

mot  useo 

tN*A* 

LIST? 

«*CSJ  T«JUTT  WE  fvsifclf  *tONC  **<  YwCWSt  *»i»  ILJMI 

*rtt^ 

LIST? 

tnoottlc  srtt«M» 

f»4V 

LIST? 

MOT  U«0 

LIST? 

n*e  «seei 

V? 

Ll?T* 

St*  Le*TL  JT*NQMiC  0*»  Tt«PsO*»U0t  1040 

V 

iMO* 

NO!  u**0 

17*4** 

l«<s* 

WOT  uLE  r 

vv 

LlS»> 

**«.L  50*1  ;  «rT?»t^i 

V** 

list? 

T*ur  *iftsreto  irT.sno 

iMC* 

V»r  VITJNStO  ITT  %t^Ut  SOI 

***«*« 

Ll«T» 

wfT  uses 

v*» 

list? 

wNruSriwe)  ccr*KiLffM*  -  ?Lur  *|B*#tfO 

If 

L©TC» 

wCT  utfo 

Ltt®» 

mst  uses 

«K.s 

lists 

*i**s  i«i?tNsti 

*■? 

LIST# 

«*«5^  »l  lUll 

»A 

list* 

*  MHv****  TO  Ll<«CljT  t'«  Ltffl* 

l«T" 

«5»  VJEO 

** 

LISTi? 

»  8SLt*w?|  tM  (*jh,  in;  till 

LtSTl* 

***S  0*  E*  l|  *!!■>  tlMf  ,11/St?; 

** 

LSS»? 

WOT  USES 

*• 

Lis** 

*  «*W«Nt  *©  TUTeCL-TIwE  Lt>S»r 

n 

LI  *.?■*■ 

?  T£-  ijN-,yT  IMt  LOO* 

?*« 

VLKT"* 

«*?*  LLT''S 

?♦ 

L|S*I* 

t  SIS**!^*  1*.  t*»I.  **|*  S*ST£«  ePT, 

L  T  STlS 

wit  o*  c**jmlu  e»  «  *i»-  ♦»«*  iff/u(i 
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SECTION  4 


INPUT  VARIABLES  AND  SAMPLE  CASE 


4.1  INPUT  VARIABLES 

Input  data  for  the  MILSTOL  program  is  read  by  the  executive  routine,  the  trajectory 
subroutines,  the  physical  data  subroutines,  and  by  the  page  heading  subroutine  (SKIP). 
Except  for  the  page  heading  subroutine,  all  inputs  are  in  namelist  format.  The 
SKIP  subroutine  reads  two  80 -character  title  cards  and  prints  them  at  the  .op  of 
each  output  page.  The  majority  of  the  variables  in  t  TAKEOFF  and  LANDING  sub¬ 
routines  are  input  via  a  data  statement  and  correspond  to  the  ground  rules  specified 
in  Figures  1-1  md  1-2.  If  these  values  are  to  be  changed,  the  data  statement  may 
be  overridden  by  inputting  the  appropriate  variable  in  the  namelist. 


Tables  4-1  through  4-6  describe  all  MILSTOL  input  variables.  The  namelist  MAIN  is 
read  by  the  executive  routine  MILSTOL,  anti  its  variables  are  described  in  Table  4-1, 
Table  -‘-2  is  the  description  of  the  title  cards  for  the  SKIP  subroutine.  Inputs  for  til? 
maneuver  driving  subroutines,  TAKEOFF  and  LANDING,  are  contained  in  the  name- 
lists  TAKEOH  and  LANDi,  Tables  4-3  and  4-4  respectively.  AEItTl  and  THT1  are 
the  namelists  for  the  aemlyuamic  aad  propulsion  data  subroutines  and  are  described 
in  Tables  4-5  and  i-6. 


4.2  SAMPLE  CASE 

The  sample  input  data  presented  in  Tabic  4-f  is  a  representative  externally  blown 
flap  configuration  with  an  aspect  ratio  of  §  and  a  quarter  chord  sweep  angle  of 
25  degrees.  The  aerodynamic  and  propulsion  data  is  consistent  with  that  used  in 
Reference  7  for  the  configuration  definition  studies.  The  input  data  of  Table  4-7  is 
set  up  to  calculate  the  ’balanced''  takeoff  d,  stance  for  a  135, 000-pound  aircraft  at 
a  runway  altitude  of’  S.500  feci  on  a  MIL-STD-210A  hot  day.  After  completion  of 
the  takeoff  calculation,  the  program  wilt  read  another  set  of  inputs  and  calculate 
a  landing  at  the  same  runway  and  weight  eouditiotis.  The  resulting  takeoff  cal¬ 
culation  output  for  this  sample  case  Is  shown  in  Table  4-@;  the  landing  calculation 
output  Is  shown  In  Table  4-9. 
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Table  4-1.  Definition  of  MILSTOL  Input  Variables 


DEFINITION  OF  VARIABLES  IN  NAMELIST  »KAIM* 

NWTa  THE  NUMBER  OF  WEIGHTS  IN  THE  WEIGHT  MATRIX 
WTLa  THE  LIST  OF  WIGHT*  <LB».) 

NHFa  THE  NUMBER  OF  ALTITUDES  IN  THE  ALTITUDE  MATRIX 
HFLa  THE  list  OF  RUNWAY  ALTITUDES  (FT.) 

VWKa  THE  weADWINO  COMPONENT  IKEA*) 

|  ATM  I*  TYC  INDEX  FOR  SELECT  I  NO  ATMOSPHERIC  PROPERTIES 
I  ATM.  0  STANDARD  CAY  TEMPERATURES  (COMPATIBLE  WITH  US  STAX?  AM) 

AND  ICAO  TEMPERATURES) 

•  |  MIL-STO-EI0A  TROPIC  TEMPERATURES 
■  l  MIL-STD'EtOA  POLAR  TEMPERATURES 
a  3  MIL-STD-210A  HOT  TEMPERATURES 
a  4  M1L-STO-2IOA  COLD  TEMPERATURES 
a  9  TEMPERATURE  IS  SPECIFIED  AS  *TEKF*  IOE0  FI 
INSEO  IS  THE  SEOUCNCE  FOR  DATA  INPUT 
iNSEOa  2  SUB ROUTINE  SKIP  -  TWO  TITLE  CAROS 

a  3  SUSROUTINE  AERO  I  -  NAHELIST  AERTl  o 

a  4  SUBROUTINE  THSTJ  -  NAHELIST  THT1 

a  9  SUBROUTINE  TAKEOFF  -  NAMELIST  TAKE  OF) 

a  *  SUBROUTINE  LANDING  -  NAMELIST  LAAOl 

a  T  TERMINATES  READING  INPUTS  AND  BEGINS  EXECUTION 

a  B  CTOP  -  END  CF  JOB 

evsco  is  the  execution  sequence  ithe  numbering  seoucnce  IS  the 

SAME  AS  INSEO I 

EXAMPLES  -  |NSE0*2.5,3,4,?,  CAUSES  INPUTS  J  BE  READ  BY 

SUBROUTINES  SKIP. TAKEOFF* AEROI.ThSTI  IN  ORDER  AND  Then 

execute  the  job 

CXSKCaS.e.  EXECUTES  SUSROUTINE  TAKEOFF  AND  THEM  CALLS 
STOP 


Table  4-2.  Definition  of  SKIP  Title  Cards. 


T*  INPUT  CONSISTS  OF  TWO  ALPHANUMERIC  TITLE  CAROS 
CARD  I  -  COLUMNS  C  THROUGH  <0  ARE  RESERVED  FOR  ENGINE 
IDENTIFICATION 

COLUMNS  t|  THROUGH  SO  ARE  PRIHTEO  AS  A  TITLE  LINE 
CARO  *  -  COLUMNS  |  THROUGH  BO  ARE  PRINTED  AS  A  SCCOAO 
TITLE  LINE 


Table  4-3.  Definition  of  TAKEOFF  Input  Variables. 


OFFIHtTION  OP  VARIABLES  IN  NAMELIST  (TAKEOFl* 

VCS  *  HAT|0  OF  LIFTOFF  SPEED  TO  AIB  N|N|NUM  CONTROL  SPEED 
OGLOa INCREMENTAL  LOAO  FACTOR  REQUIRED  AT  LIFTOFF 
DFLPaFLAP  SETTING  IN  DEGREES 
RCOEF  a  ROLLING  COEFFICIENT  OF  FRICTION. 

BCOEF  a  COEFFICIENT  OF  BRAKING  FRICTION 

VMCGKa  MINIMUM  CONTROL  SPEED  ON  THE  GRODNO  (KNOTS! 

VMCAK  a  MINIMUM  CONTROL  SPEED  IN  THE  AIR  IKNDTSI 
POTATN  ROTATION  RATE  (DEGREES  PER  SECOND) 

T|MF  a  REACTION  T|HR  FOR  ENGINE  FAILURE  ISEC) 

T|PB  «  BRAKING  DELAY  AFTER  *TIMR*  (SEC) 

OOTPT  EQUAL  TO  ZERO  SUPRESSCS  PRINTING  OUTPUT  FOR 
SEGMENTS  A, 6  AND  s. 

THE  FOLLOWING  VALUES  ARE  ENTERED  AT  TIME  OF  LOAO  INC  ANO  ARE  USED 
UNTIL  OVERRIOEN  BV  READING  THE  APPROPRIATE  VARIABLES  IN  TAKEOFl 
VSCal.OS 
DGLOaO.IO 
CFLR«F9, 

RCOEFaO.IO 

BCOCF.O.JO 

TIMRaJ.O 

TIMB«2.0 

ROTATNaB.O 

ROTRTa), 
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Table  4-4.  Definition  of  LANDING  Input  Variables 


0FFINIT1CN  OF  VARIABLES  IN  NAMELIST  •  LANS  I  » 

oelp»  elar  setting  ioegi 

ME  ON.  OBSTACLE  HEIGHT  ITT  ) 

ROTATNs  ROTATION  RATE  (DCG/SECI 

VMCGKa  MINIMUM  CONTROL  SPEED  ON  THE  GROUND  IXTSl 

VMCA<.  MINIMUM  CONTROL  speed  IM  TIC  AIR  «KTS> 

TERM.  TIME  delay  AFTER  TOUCN00VN  eor  brake  application  tseci 

TSP*  TIME  delay  AETER  TOUChOOWM  eor  spoiler  deployment  tstci 

TaFVf  TIME  DELAY  AETr.R  TOUCHDOWN  EOR  THRUST  REVERSAL  ISECl 

APR*  RATIO  OE  APPROACH  SPEED  TO  MINIMUM  CONTROL  SPEEO 

OGTlia  INCREMENTAL  LOAO  EACTOR  AVAILABLE  AT  THE  OBSTACLE 

OGTOa  INCREMENTAL  LOAO  EACTOR  AVAILABLE  AT  TOUCHDOWN 

RCOEE.  ROLLING  COEFFICIENT  OE  FRICTION 

BCOEF.  BRAKING  COEFFICIENT  OE  FRICTION 

gamma.  initial  flight  path  angle  ioegi  (negative  is  oescenoingi 

R«.  MAXIMUM  RATE  OE  SINK  AT  TOOCHOOWN  I FT/SECS 

I POSITIVE  IS  OESCENOINGI 

ROTPT*  EQUAL  TO  ZERO  SUPRESSES  PRINTING  OF  THE  ENTIRE  TIME  HISTORY 

The  FOLLOWING  VARIABLES  Are  ENTERED  AT  TIME  OE  LOADING  AND  ARE  USED 
UNTIL  OVERRIDDEN  BY  READING  THE  APPROPRIATE  VARIABLES  IN  LAND1 
OELP  ■  60.0 
ME OB  *  SO.O 
rotaTn  .  e.o 
TSRK  .  2.0 
TSP  ■  2.0 
TREV  «  2.0 
APR  *  !*J0 
DOTH  *  0*30 
OGTD  »  0.1* 

OCOFF  *  0*|0 
6C0CF  •  0*30 
GAMMA  *  -?*50 
«S  •  10,0 
ROTPT  «  WO 


Table  4-5.  Definition  of  AER01  Input  Variables. 


DEFINITION  OE  VARIABLES  IN  NAMELIST  *AERTI  t 
S  ■  WING  REFERENCE  AREA, 

BW  a  WING  SPAN 

CAHT.  HEIGHT  OE  THE  QUARTER  CHORD  ABOVE  THE  GROUND  irT| 

GALMW.  ANGLE  OE  ATTACK  EOR  GROUND  CONTACT  (OEG| 

alpmx  •  maximum  allowable  ancle  of  attack  -  usually  a  cl  limit 

(OEG| 

A CL MO  a  STALL  ANGLE  OE  ATTACK  (DEG! 

EPP  a  EQUIVALENT  FUAT  PLATE  AREA  OE  ADDITIONAL  ORAG  ITEMS  ISO  ETl 

winco  a  vine  incidence ioegi 

OCLSP  a  INCREMENTAL  LIFT  COEFFICIENT  DUE  TO  LIFT  DUMPERS 

ccnsp  a  incremental  drag  coefficient  due  to  lift  dumpers 

EOOCL  a  INCREMENTAL  LIFT  COEFFICIENT  DUE  TO  ENG  OUT  CONTROLS 

foocd  *  incremental  orag  coefficient  due  to  eng  out  controls 

NX.NV.NZ  APE  THE  NUMBER  OF  XC.VC.ZC 
XC  a  THE  LIST  OF  ANGLES  OF  ATTACK  IOEGI 
VC  a  THE  LIST  OF  THRUST  COEFFICIENTS 
ZC  a  THE  LIST  OF  flap  DEFLECTIONS  I DEG I 

CLA  a  lift  COEFFICIENT  AS  A  FUNCTION  OF  I ALPHA .CT.FLAP  DEFLECTION! 
COA  a  ORAG  COEFFICIENT  AS  A  FUNCTION  OF  (ALPHA. CT.FLAP  DEFLECTION! 
.TIG  a  I  FOR  MECHANICAL  FLAPS  (X.U5  VECTORED  THRUST 

2  FOR  EXTERNALLY  BLOWN  FLAPS 

3  FOR  INTERNALLY  BLOWN  FLAPS 


uuuou 


Table  4-6.  Definition  of  THST1  Input  Variables, 


e  definition  of  var iables  in  namelist  »thts • 

C  H-  NUMBER  OF  X.V  POINTS  IN  CALM  TABLE 

c  ws-  the  velocity  table  r«  the  propulsion  tables  in  etas 
e  tuseo  as  ?>c  independent  variable  in  all  tables) 

C  VHST-GROSS  THRUST  TABLE  AT  MAX  MOVER  I  IN  LBS.) 

C  HORO-RAM  CM  AO  TABLE  AT  MAX  MOVER  UN  LBS.) 

C  TICL-6R0SS  THRUST  TABLE  AT  IDLE  MOWER  (IN  LBS.) 

C  TSLT— GROSS  THRUST  AT  T«  SLOT  EXIT  1 1ST  CONFIG! I  AT  MAX  MOVER  (LBS) 

C  TRrV-MAX  REVERSE  THRUST  (EXPRESSED  AS  A  NEGATIVE  VALUE  -  IN  U)S.> 

C  OVNF-W I  SCHILLING  DRAG  FOR  A  OSAO  ENGINE  I  IN  LBS.) 

e 

C  THiO-THRUST  VECTOR  INCIDENCE  REM.  TO  A  HATER  LINE  IN  DEG. 

C  ENCNO— THE  NUMBER  OF  ENGINES 

C  SCALE-SCALING  FACTOR  FOR  THE  PROPULSION  DATA 
C  NTNO»0  NO  REVERSE  TtCUST 

C  *1  ALL  ENGINES  REVERSING 

C  «2  ENGINE  OUT  REVERSING  PROCEDURE 

TK  FOLLOWING  VARIABLES  ARE  ENTERED  AT  TIME  OF  LOADING  AND  ARC  USED 
until  OVERRIDDEN  8V  HEADING  THE  APPROPRIATE  VARIABLES  IN  TNT I 
ENCNO  •  4.0 
SCALE  •  1.0 
NENG  a  | 


Table  4-7.  Sample  Input  Data. 


AMAIN  NVTal.  VTL« 133000,  ,  NNF*t  *  H?L»SS0O.«  I  ATM.  3,  VHK.O.. 
INSC0a2,S,3,4,7,  EXSC0.9.T.  • 

H  I3-F2B  ENGINE  SAMPLE  CASE 

FOR  AN  EXTERNALLY  BLOWN  FLAP  CONFIGURATION 
•TAKEOFt  VSCM.OS,  0FLP*2S,»  RCOEF.O.I,  BCOEF.O.3.  VMCGK.6B. >  VMCAJC-BB,. 

POTATN.fi,  •  TIMR.)..  TIMB.i..  DCLO.O.U  ROTPT.O..  S 
6ACPTI  S.lSSOt,  ALPMX.20, «  GALMX.20, •  AClMO«ZT,,  JFIC»*»  DCLSMo-2.0, 
OC0SP.0.2S,  E00CL*-0,43,  eooco.o.ob* 

Nit,  I,,  i  r^Y*4,t  I2.A .  ; 

XC*  —A, ,0, ,4, ,8. ,12, *16,  ,20»»24, •2B,,32«, 

YC«  0,0. 1.0. 3.0, 3.0. 

2Ca  0,0 .30.0 .40 .0,60.0, 

CL  A. 


-0.29, 

0.03, 

0,35, 

0,67, 

0.99. 

1,30, 

1.42, 

1.9A, 

2.28, 

2*90, 

-0.33. 

0.01, 

0,34, 

0,74, 

1  *)A, 

I. 31. 

1,91. 

2.30, 

2,70, 

3.10, 

-0.73, 

-0.25. 

0,33, 

0,86, 

1.39, 

1.93, 

2. AT. 

2,99, 

3.82, 

4.00, 

-1.00, 

-0.33, 

0,31, 

0,96. 

1.63, 

2,3t« 

2,04, 

3. 64, 

A, 3|, 

4,90* 

O.BS, 

1.19* 

1.31, 

1,83, 

2*19, 

2.32, 

2,82, 

3.09, 

3,  IS, 

3.00, 

0.03, 

l  ,AA, 

1.96, 

2.49, 

3.01, 

3.33. 

A. 03, 

A. AS, 

4.63, 

4,80, 

1.30, 

1.98. 

2,64, 

3,33, 

4,02, 

4,72  , 

3.39, 

6,07, 

6,71. 

7,10, 

1.30, 

2.33, 

3,  |4  , 

3.93, 

4,73, 

3.53, 

6.28. 

7,11, 

7,67. 

B«90, 

1.33, 

1,63, 

1.96, 

2.28, 

2.89, 

2.90, 

3.21* 

3.30, 

3.33, 

3,00, 

2.30, 

2.9)  , 

3.37, 

3,81  , 

4,29, 

4.74, 

3,18, 

8,45, 

8.26, 

8.00. 

3.80, 

4, A3, 

8,03, 

8,66, 

6,28, 

6.63. 

7,30, 

7,68, 

7,90, 

§.00, 

*.ao, 

3,43, 

6.12, 

6, BO, 

7.40, 

4.13, 

6.8T, 

9.11, 

9,42, 

9.80, 

1.38, 

1  .AS* 

1,96, 

2,2B, 

2,39, 

2.40, 

3,21. 

3,80, 

3.33, 

3,00, 

2.30, 

2.91* 

3.37, 

3. at, 

4,29, 

A.  74, 

8. IS, 

8,48, 

8,26, 

8,00, 

3, 80, 

A, A3, 

8,03, 

8,64, 

6.20, 

4.83, 

7,30. 

7,68, 

7,90. 

a, co, 

A.ao, 

3,  A3, 

6.12, 

6, BO, 

7.48. 

S.I3* 

4.87* 

9.11, 

9.42, 

9.80, 

^7v.:ra  n  TvrS-srt-i/'?^’}  >  ?rV?V;!  7.  ’.li.vf  t 


ww*wep*«  •■wow  ;•» 


Table  4-7.  Sample  Input  Data.  (Contd) 


0.1S, 

0,19. 

0.18, 

0.16, 

0.19, 

^  ; 

0.20, 

0,22, 

0.26, 

0,33. 

0,49, 

f  \ 

-0.90. 

-0.90, 

-0,88, 

—0.85, 

-0,82, 

-0.78, 

-0.73. 

-0,67, 

-0,98, 

— o,*0. 

•(. 

-2.90, 

-2,90, 

-2,89, 

-2,63, 

-2.79. 

-2.73, 

-2.66, 

-2,55, 

-2.41,' 

-2.10. 

k 

-4.90, 

-4.90. 

-4.89, 

•4,83, 

-4.77, 

A 

—*•68, 

-4.56, 

—4,43, 

-4.23, 

-3.90. 

0.21, 

0.21, 

0.24, 

0.30, 

0.38, 

0.47, 

0.98, 

0,7t « 

0.91 . 

0.69, 

7 

•0,4*4 

-0,42, 

-0,36, 

•0,26,. 

-0.12, 

0.08, 

0,32, 

0.94, 

0.83, 

0.99, 
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PROGRAM  AND  SUBROUTINE  SOURCE  LISTINGS 


The  following  source  listings  are  contained  in  this  appendix. 
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Description 

Page 

MILSTOL 

Main  Program 

1-1 

MILSTOL 

Overlay  Dump  Program 

1-3 

AEQFM 

Equations  of  Motion  Subroutine 

1-4 

AEROl 

Aerodynamic  Data  Subroutine 

1-7 

ALIN 

Linear  Equation  Function 

1-9 

ATMOS 

Atmospheric  Properties  Subroutine 

1-10 

FIND 

1-Dimensional  Table  Lookup  Subroutine 

1-12 

GILL 

Integration  Subroutine 

1-13 

GINTG 

Integration  Driver  Subroutine 

1-14 

HEAD 

Page  Heading  Subroutine 

1-15 

INTP2 

Curve  Fitting  Subroutine 

1-15 

KABD 

Hyperbolic  Curve  Fit  Solution  Subroutine 

we 

LAGRA 

Lagrantan  Interpolation  Subroutine 

i-ia 

LANDING 

Landing  Trajectory'  Driver  Subroutine 

1-17 

LOOK 

3-  Dimensional  Table  Lookup  Subroutine 

1-25 

NWRP2 

Newton- Wrapsoc  iteration  Subroutine 

1-25 

SKIP 

Page  Eject  Subroutine 

S— 30 

SL 

Linear  Slope  Function 

1-31 

SMLT2 

Simultaneous  Equation  Solution  Subroutine 

1-32 

TAKEOFF 

Takeoff  Trajectory  Driver  Subroutine 

1-33 

THST1 

Propulsion  Data  Subroutine 

1-41 

HUM 

Aircraft  Trimming  Subroutine 

1-43 
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PROGRAM  MIL,ST0u« input .OUTPUT, TAPES* INPUT »TAPE6»0UTPUT ) 

C  COC  6400  SHORT  TAKEOFF  AND  LANOING  COMPUTER  PROGRAM  USING  MILITARY 
C  GROUND  RULES 

C  THIS  PROGRAM  REQUIRES  INPUTS 

C  THIS  PROGRAM  CALLS  THE  FOLLOWING  ENTRIES 

C  SK1P.TAKE0FF.LAND«NG.AER0I«TMSTJ 

INTEGER  CARD* PRINT .PAGE  »EXSEG 
PFAL  LIFT. MACH, NU.KTOF.NMTOF 


OIMFNSION  EXSFOI2), INSE0I3I.WTLI 1 0 ) ,HFL<10> 

COMMON  /LIST1/LINE. PAGE. LIMIT. CARD. PRINT, INP, I PUNCH. I DATE, HD 160 > 
COMMON  /L I ST2/0EL YD , ALFAR , THRTL  » GAMR , M 1 ND  t  VTF • WT , HF , OT , T I M$ , XF , 

*  FUEL , MACK, VWF  «L I FT .DRAG*  TKRST . RCF .ENQNO, THV *  I N 

COMMON  A.  I  ST  3/CL  »  CO  «  S  ,  Q5  »  TH  |  R  .  /.L  THR  .  CLMAX  ,  AR  ,  CLAR  ,  CMACG  «  ALPMO  ,  CX  , 

*  C2,0,IMOM,WINCP.CMDMP.OSC,ALMXR,VSF,8W 
COMMON  /L I ST4/FF5 , UPQW , TGROS . DW A « THMOM , TMREO 

COMMON  /L I  STS/S I G  »  SOUND , Nu , T  EMR  «  P  AKB . I A  TM , TEMF , DSODH , ORHO 

COMMON  /L I ST6/RT0D ,DTOR ,KTOF , FTOK  »  NMTOF • FTONM , RHOZ , RH022 , GZ ,P2 . TZ 

COMMON/L I ST 8/ I NOEX «COEF«GN. GT , AN. AT , AX, AZ , I TRM, VUPPEW 

COMMON/L  f  ST9/ANS , ANS2 , NO , XA . YA , ZA . I ND I C . 1 FL AG ,NOER 1 3  > , OBOX , D3DY  * 

*  PP02 ,DB20X ,0920Y  <  082DZ 

COMMONA  I  ST  I S/  r  ALMX,  I  NT  ,  X£DOT  ,  ZEDOT  .XE  .20.  ILOOP  »DT  |  ME  ,  I  AT  ,F'_  |  FT  , 
* OFLP  «  C4HT .GALMXR 

COMMON/L GEOM/SMSW , SCR , QQOT .RTH.OTnTM.RWA.OTmwa.OSThZ »  0 WLH2  *  I YY . 

1  XCG«ZCG,C0AR.PFN,SINTM,COSTM,ANGLE,S3NAL,COSAL,RHZ,OTmmZ,VTSQ, 

2  L "HT .MOOT ,ThETR,U.W 
C0m.'.=om/C0NTR0L/JF|G»  IREV,  ISP.NENG 
COMMON  a.  I ST99/ 1 ERR 

0 AT  ART0N , KTOF , NMTOF , RHOZ ♦ RH022  »  GZ ♦ »Z . T  2 /5  7 . 29S 7795 . I ,6uy . 60  76 . b , , 0 
lO?3^68R«u4,,oei  1604462. 32*  I74C*9,2j  16.2218,510,67/ 

DEFINITION  O-'  VARIABLES  IN  NAMELIST  «MA|N* 

NW'T»  THE  NUM6“»  OF  WEIGHTS  IN  THE  WEIGHT  MATRIX 
WTL*  THE  LIST  OF  WEIGHTS  ILSS.l 

NMF»  ThE  NUMQER  OF  ALTITUDES  IN  TmE  ALTITUDE  MATRIX 
MFL»  THE  LIST  OF  ..UNWAY  ALTITUDES  IFT.J 
WK«  TmE  HEADWIND  COMPONENT  (KSASl 

14TM  IS  THE  index  for  SELECTING  ATMOSPHERIC  PROPERTIES 
lATM*  0  STANDARD  pay  TEMPERATURES  tC.OMPAT  IDLE  WITH  US  STANDARD 
AND  ICAO  TEMPERA TURESI 

•  l  MIL-ST0-2IGA  TROPIC  TEMPER/ TyRES 

•  2  MJL->*TO~2tOA  POLAR  TEMOFRATutfES 
w  3  MIL-STD-2I0A  «0T  TEMPERATURES 

«  *  m|L-STO<-25FA  CCSLQ  TgMPeRATLWEC 

4  9  TEMPERATURE  IS  SPEC {FIFO  AS  » Tfwr »  COEG  f> 

iNT.ro  is  tmf  sequence  fd»  data  input 
tN4t0*  2  SUBROUTINE  SKIP  -  Two  TltLS  CAPOS 

4  3  SUBROUTINE  escoi  -  nakelIST  aCRT i 

4  *  SUBROUT  I Mg  TMStt  -  NAMELIST  TnTI 

4  9  SUDROUTINS  TAxEOFF  -  NAMElIST  Taweufi 

•  6  SUBROUTINE  landing  «  NA^LSST  land. 

•  T  TERMINATES  READING  |«A/T j»  AND  BEGINS  EXECUTION 
4  «  STOP  -  E*D  OF  UC0 


o  u  o  u  u  u  u 


FX<FO  IS  THF  EXECUTION  SEOUFNCF  (THE  NUMSERJNG  SEQUENCE  IS  THE 
SAME  AS  INSrO) 

FXAMPi.es  -  !NSE0*2.5.3.4,7,  CAUSES  INPUTS  TO  BE  BEAD  BY 

SUBROUTINES  SKfP, TAKEOFF. AEOQl.TMST!  IN  ORDER  AND  THEN 
EXECUTE  THE  JOB 

EXSE0*5«8<  EXECUTES  SUBROUTINE  TAKEOFF  AND  THEN  CALLS 
STOP 

NAMEL I  St/MA  IN/NWT «  WTL«NHF .  fs? L  ♦  VWK .  |  ATM « T£MF «  I NS£0«£X5EQ 

C/s90»5 
POINT  «6 

ir*ATe«MioaF»3i9) 

L|M|T*55 
PTO»*I*/»TOO 
FTOK*  I  ,/K.TOF 
FTONM* I ./NMTOF 
PAGE  *  I 
I  IN*  I 

OTAOCCARO.MAIN) 

VPITEIPfilNT.MAINJ 
f  •! 

5  iNO.jNSEOm 

GO  TQl99«20i30*4''«F0*60i  ISO*  |6D>  •  |NP 
20  CALL  SKIP 
GO  TO  90 
30  CALL  AE«0| 

GO  TO  90 
AO  CALL  TMSTl 
GO  TO  90 
00  CALL  TAKEOFF 
GO  TO  93 
60  CALL  LANDING 
60  TO  90 

99  woiteippint.ioooi 

IOOO  rpevATt  {Ml  •  lex.AC-Meatfoo  IN  NAMELIST  *ha|n»,  INSEO  OR  £xs€Oi 

Kfpp 

OO  |F(|N.GT,| »GO  to  I AO 
!•!♦» 
f-'  TO  5 
150  »»«S 
IN  •  Z 

»*J  IF**e*Sfc'OM  n~* 

PO 

PO  I S3  *®|»NST 
Jl 

CALL  ATMOS 

V**'*VW«TQA/L89TIS|6l 

WT«%fTL»K| 

im«!> 

i*jt»*e.xsF9«  i  n 

GO  TO  iSO.GO* ! . $6S> . IF* 

I  40  CONTINUE 
ISS  rpNTfNvC 
IS*  CONTINUE 

t  **n*s 

go  to  m 

1*9  vPttf  rP3|NT«t*?l ) 

»G9I  )H|.50K,|CHfM5  OF  JC#  J 


C  »  M  t>  U  N 


»M*:.^?*!;-*>J^>  V  ‘"jf  *^0  JSf^ViJ^V'^SFK  ;_•_'  >  :“• 


PROGRAM  MIlSTOLI l NPUT, OUTPUT .TAPES* INPUT, TAPE6«0UTPUT> 

C  CDC  6400  DUMP  CVFBLAV  PROGRAM  POP  »M!LSTCL» 

C  THP  OUMP  IS  A  FLOATING  POINT  NAMELIST  PRINTOUT 

INTEGER  CAHD»PR|NT,PAG£,EXS£0 
REAL  LIFT. MACH, NU.KTOF.NMTOF 

COMMON  /LIST! /LINE, PAGE.LIMI T, CARO, PRINT, |NP, IPUNCH, 10ATE,HO(60) 
COMMON  /L I ST2/0CLT0 , ALF AR , TmRTL .GAMR,W|NO»VTF,WT.HF.OT,TtMS*XF, 

*  FUEL.MACH.VWF, LIFT, DRAG, TMRST.RCF.ENCNO.TMV. IN 

COMMON  /LIST3/CL«CO,S,OSiTMlR,ALTHR,CLMAX,AR,CLAR,CMACG,ALPHO,CX, 

#  C2 ,0 , I MOM , V I NCR , CMOMP , OSC , ALMXR  »  VSF • BW 
COMMON  /L I ST4/FFS, JPOW , TGB0S  *  OX A , THMOM , THREO 

COMMON  /L I ST9/S I G , SOUND  «Nu , TEMP , P AMB » I ATM , TEMP . DSOOM , ORHO 

COMMON  /L I ST6/RT0C ,OTOR,KTOF .  FTOK , NMTOF , FTONM , RM02 « RW022 » GZ  »PZ  ,  TZ 

COMMON/L I  STB/ 1 NOEX , CO€F , GN , GT , AN, AT , AX , AZ , I TRM , VUPPER 

COMMON/L I ST9/ANS , ANS2  *  NO ,  XA  ,  YA  ,  ZA  , INOtC, IFLAG.NQERO) *Q60X,OBOY* 

•  OBOZ , 0820X , 0920 Y • 0B2DZ 

COMMON/L I  ST 1 5/ I ALMX . INT . X£OOT . ZEOOT , XE , 2£ , ILOOP ,0T I M£ ,  l AT ,FL I FT , 
1 OFL  P»C4MT,GALMXR 

C*  •  *ON/LGEOM/ShSW»SCR*OOOT,RTH,OTmTm,R«A,OTm*!A»DSTh2,DWLm2.  IVY, 

1  x  C.ZCG.CBAR.PFN, SI NTm,C05Tm» ANGLE »S»NAl.CGSAL.Rh2.QThkZ.VTSG. 

2  UOOT, MOOT, the tr,u,w 
COMMON/CONTBOL/JFIG, {REV. 1SP,N£NG 
COMMONA.  I ST99/ 1  ERR 

NAMEL I ST/OUMP/l  | N£  .PAGE  . L  1 M I  T  .CARO , PR t NT  ,  I NP ,  | PvNCM ,OFLP , 

1  OELTO,ALFA«.THRTL.GAMR,MINO,VTF,»fT.MF.OT.TJMS,XF. 

Fuel  ,  HACM  ,  V*F  «  L I F  T  ,  OR  AG  ,  TKRST  .  RCF  ,  ENGNO  ,  THV  ,  I N  .  C  «k  HT  •  G  ALMXR 

»CL«CO»S,qS»TM|R»ALTMR.CLMAX,AS,CLAR,CMACG»ALPHO,CX, 
C2.0.,»O-..,~rS.C»KP.^l.U-»=.^5a.l>(ii£0 

.  S I S  ,  SOWO  ,  NU  .  TEMO  ,  P A.M0  ,  I A  TM  ,  TE WF  ,  DSOOM  ,  ORHO 
, RTOO »OTCR  .KTCF  .FT OK  , NMTQF  .FTONM, RmOZ  , RM>2*» , G2  ,PZ  .  T z 
*  INOEX.CCEF .GN.GT «  AN,  AT,  AX ,  AZ •  I  TR*», vuPPEft 
,AN^,ANt?.2,N0,XA.VA.ZA.  INOtC*  ;.-LAGiNOER  .DS'DX.OODY, 

I  OBOZ . 0B29X .CK»20V , 0R20Z , I AL«X .INT. x£OQT . ZEOOT . X£ . Zt . 1 1  OOP .0  T J  Mg , 

2 1 At ,FL l FT  «  SMS* . SCR , GOOT «RTM , DYhTm ,R«A . CTh*  A • OS  TMZ . D*L.hJ , I vy , 

3  XCG.ZCG.CB-S.PFN.St^TM.CCSTM. ANGLE »S!NAl,CCS*L.»hZ,OTmm2.VTSO. 

*  UOOT.xOOT.THFTC.U.Kf.OFtG.  IBEV.  ISP.  IERO 

%OITEIB.0UMP» 

STOP 
FNP 
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subroutine  aeofm  J 

C  6400  SUBROUTINE  FOR  CALCULATING  ACCELERATIONS.  J 

J 

c  this  subroutine  calls  for  the  following  entries.  j 

C  THSTI.  AERQ2 .  TRIM  J 

J 

C  ALL  index  OPTIONS  ASSUME  WEIGHT  AND  ALTITUOE  ARE  SPECIFIED)  J 

C  ALL  EXCEPT  3  ASSUME  GAMR  SPECIFIED.  J 

C  ALL  EXCEPT  6  ASSUME  POWER  SETTING  SPECIFIED.  J 

C  ALL  EXCEPT  7  ASSUME  VELOCITY  SPECIFIED.  J 

C  6N  AND  6T  ARE  RESPECTIVELY  THE  NORMAL  ANO  TANGENTIAL  LOAD  FACTORS  J 

C  INCiEX  a  j  FOR  CALCULATING  GROUND  ACCELERATIONS! GIVEN  ALFAR.  GAMR.  J 

C  AND  CO£F> .  j 

C  INDEX  *  2  FOR  CALCULATING  AT  AND  AN! GIVEN  ALFAR.  GAMR).  j 

C  INDEX  ■  3  FOR  CALCULATING  ALFAR  AND  6AMR t G I V£N  GN  AND  GT).  J 

c  Index  *  «  for  calculating  alfar  and  anigiven  gamr  and  gt>.  j 

C  INDEX  »  5  FOR  CALCULATING  ALFAR  AND  ATiGIVEN  GAMR  AND  GN).  J 

C  INDEX  a  6  FOR  CALCULATING  ALFAR. ThRTL*  IGIVEN  GAMR.GN.GT).  J 

C  INDEX  a  7  FOR  CALCULATING  ALFAR «VTF (GIVEN  GAMR.GN.GT).  J 

J 

INTEGER  CARO. PRINT, PAGE  j 

real  LIFT.MACM.NU.KTCF.NNTOF.IYY  j 

DIMENSION  SA!3).SV!6).SPV(G) ,ST(3) ,UA<3) . JV(3) .A! 6) .V(5) »PV16) .T!6  J 

I)  J 

COMMON  /LISTI/  L 1NE.PA6E  *L  1M  J  t  .CARO.  PR  I  NT  »  I  NP,  I  PUNCH,  |  DATE  .HO!  60 1  J 

COMMON  /LSST2/  OELTO. ALFAR, THRTL.GAMR. MIND. vTP.wT.MF.OTtTIMS.XF.FU  J 
IFL.MACH.VWF.LIFT.ORAG.TmRST.RCF.ENGNO.TMV.IN  J 

COMMON  /L I  ST  3/  CL.CD.S.OS.TmIR.ALThR.CLMAX.AR.CLAR.CMACC .ALPMO.CX.  j 
IC?.Ot  tMOM.WINCR.CMQMP.QSC.ALMXR.VSF.aw  J 

COMMON  /l l ST4/  FFS.jROW.TGROS.OWA.TmmOM.TMREQ  j 

COMMON  /LIST'S/  s  I G.  SOUND.  NU.TEMW,  PAM0,  i  ATM,  TEmF.dsoDm.ORmQ  j 

COMMON  A.IST6/  RTOD.DTOR.iCTOF,FTOK.NMTOF,FTONM,RmOZ.Rn022,G2.P2.T2  j 
COMMON  /<_  \  STS/  INOSX.COFF.GN.GT.AN.AT,  AXE.A2E,  ITHM.VUPPER  J 

COMMON  /LOfeO**/  ^MSw.SCR»0OOT.RTM.DTMTM,RWA%aTMAA.DSTM2.0wUM^,  IW.JC  j 
I  CO «?CC. « CBAH  , MFN ,  S I NTM . COSTM ,  ANGLE  ,  S I NAL  • COSAL  •  RM£ »DTmm2  .  vTSO .UOQT  •  J 
2WOOT.TMETO.U.W  j 

COMMON  /L ISTQO/  Ifca  J 

EQUIVALENCE  tSAt | ) .OCDOT) ,  !$A!2),ACN).  CSAi3).ACTi*  lSvt|).C€LTO)  J 

OAT*  ST.SPV/.0CC3.2w.0CC2,.2,.0C3*RO6S^4,-*02..003A9OPSSA,.2»e./  J 

RfCPM* G?/WT  j 

eMT'S.RMQ?2asiG*S  J 

k»-I  J 

NV'TP.O  J 

f».fTF»0,ep5aLOUND  J 

MSYgPwO  J 

J 

C  5FT  SU&SCOfPTS  FOR  A.V.PV.T,  j 

J 

I»1  J 

IF  ISMOM.EO.3j  GO  tc  10 

JA  »  S I  *  t  J 
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L 

1 


10 

20 


30 


40 

50 


UV(t)«l 

IF  (ITRM.EQ.2)  JV( J ) «5 
IF  (INDEX.LT.3?  GO  TO  50 
J«!  +  l 
GO  TO  20 

IF  (IN0EX.LT.3)  GO  TO  100 
JA< ! )«2 

IF  ( INDEX. EQ.4 )  JA(I)»3 
JV( 15*2 

IF  (INDEX-4)  50,50,30 

IF  (INDEX.NE.6)  GO  TO  4C 

I-I  +  l 

JA( I )«3 

JV( I >«3 

JP0W«3 

CALL  THST1 

ThRNP*THRTL 

JP0W«5 

GO  TO  50 

IF  { I NDEX.EQ. 7 )  NVSTP«0 
NSOL-I 


STORE  VALUES  OF  INDEPENDENT  VARIABLES  (V),  PERTURBATION 
INCREMENTS  (PV)  AND  TOLERANCES  (T). 


60 


DO  60  I « 1 ,NSOL 
J«JV( I ) 

V ( I )«SV( J) 

PV( I )*SPV( J) 
J«JA( I ) 

T (  I  )«ST( J) 


START  CONVERGENCE  LOOP. 


IF  (INDEX-2)  1 00  <100, 70 
•*0  AGT-GT 

AGN*GN-1 , 

IF  (IN0EX.NE.6)  I TRM* 1 
IF  (INDEX-4)  1)0,80.90 
80  AGN*0.0 
GO  TO  110 

90  IF  (INDEX.EQ.5)  AGT«0.0 
GO  TO  110 
100  AGT«0,0 
6GN«0 , 0 

110  IF  ( I NDEX.NE .7 )  GO  TO  120 
OS«RHOS*VTF»VTF 
OSC«OS*CBAR 
MACH»VTF/SOUND 
120  althr«alfar-*incr+thir 

S INTH»S IN ( ALTNR ) 
COSTH»COS(ALTHR> 

CALL  THST1 
CALL  AER02 

TMIND«TGROS#COSTH-DRAG-DWA 
IF  (INDEX.NE.3)  GO  TO  130 
SNGAM.TM I ND/WT-AGT 
GAMR«AS I N< SNGAM ) 

GO  TO  140 

130  SNGAM«SIN(GAMR) 

140  CSGAM«COS(GAMR) 
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j 

j 

j 

j 

j 

j 

j 

j 

j 

j 

j 

j 

J 

j 

j 

j 

j 

j 

j 

j 

j 

j 

j 

j 

j 

j 

j 

j 

j 

j 

j 

J 

j 

J 

j 

j 

j 

J 

j 

j 

j 

j 

j 

j 

j 

j 

j 

j 

j 

u 

j 

j 

j 

j 

j 

j 

j 

j 

j 

j 


m 


1 

■Sf.A 


is 


f 


$s 


J 


'*5 


i9f 


£'  a  iy!,  wr/os.,.*.. 


re  J  FT  «e  l  FT-*TGROS*S  l  ntm-wt#csgam 
ANbFL I FT *PECPM 
IF  (INOEX-1)  150 «  15C  »  I7C 
j 50  IF  (FLIFT)  160. 170. 170 
160  frict«coef*flift 
00  70  180 
170  FR!CT«0«0 

180  FTHST«TMtN0-»FBlCT-WT»SN6AM 

at»fthst*recpm 

ACNfAN/GZ-AGN 

ACT«AT/GZ-AGV 

IF  { JMOM.EQ.O 1  60  TO  200 

IF  (JTRM.EO.O)  GO  TO  190 

ANGLE  *  ALF  AR-W I NCR 

ThETR-ANGLE+GAMR 

5INAL»SIN( ANGLE) 

COSAL»COS I ANGLE ) 

U«VTF*COSAL 

WnVTF*SINAL 

190  000T«AN*SINAL4-AT*C0SAL-Q*w 
WOOT« AT*S 1 NAL-AN*COSAL*Q*U 
AlPMO* (U*WOOT-W*UOOT ) /VTSO 

CMALD«“1 1 .94252*1 • 4008279-7. 52341 2/ ( ALFAR*RT0D-25 .84259) **2 ) 

A1-0.5*CBAR/VTF 

CMOMP- A l * ( CMQ»0+CMALD* ALPHO ) 

QCDOT« ( THM0M4 ( CMACG+CMDMP )*OSC ) /I /V 
200  IF  (ITRM.EO.O)  GO  TO  30C 

C  STORE  CURRENT  VALUES  OF  ACCELERATIONS  (A). 

00  210  I-l.NSOL 
J«r.lA  <  I  ) 

c  CONVERGE  values  OF <  A )  BY  VARYING  VALUES  OF  CV). 

CALL  TRlf  (NSOL.A»V.PV«T«K) 

IF  (K-NSOL)  230.230.220 
220  IF  ( INDEX. NE. 7)  290. 2S0 
C  STORE  RESET  VALUES  OF  INDEPENDENT  VARIABLES. 

230  nstep«nstep+i 

IF  INSTEP. GT. 33)  GO  TO  280 
C  RF SET  VALUES  OF  !V>. 

DO  240  I-l.NSOL 
J»JV( I  ) 

240  SV(U)«V<  t ) 

go  to  no 

250  IF  ( ASS (ACT) »LE , ST ( 3 )  >  GO  TO  290 
tF  (NVSTP.GT.O)  go  TO  270 
IF  (ACT.LT. 0.0)  GO  TO  270 
260  vtf.vtf+dvtf 

NSTEP-C 

K«-l 

GO  TO  no 

270  OVTF»“ACT*DVTF/( ACT-ASAV) 

asav-act 

nvstp-nvstp+i 

IF  (NVSTPeGT .12)  GO  TO  280 
GO  TO  260 
280  1FRR»3 
290  I TRM » 0 
300  0D0T*QCD0T 

rcf«vtf*sngam 

axe»at*csgam-an#sngam 

A7E»AN*CSGAM+AT*SNGAM 

return 


'Zil'fVuiyi  A, At.  £!\\  : 


c 

c 

c 


SUBROUTINE  AEROl 

6400  FORTRAN  IV  PROGRAM  TO  STORE  AND  RETRIEVE  AERODYNAMIC  DATA 
FOR  USE  WITH  TABULATED  AERODYNAMIC  DATA 

CL  AND  CD  «  F< ALPHA. THRUST  COEFFIC I ENT « FLAP  DEFLECTION) 


THIS  PROGRAM  CALLS  THE  FOLLOWING  SUBROUTINES* 

KABD 

LOOK 

INTEGER  CARD* PR I NT, PAGE 

REAL  LIFT.KT0F,MACH*NU.NMT0F* 1YY 

DIMENSION  XC< 15) *YC(5) «ZC<5) «CLA(375) *CDA(375) 

common  /LISTl/  L IN£«PAGE«L IMIT ♦ CARD* PR  I  NT  « INP, IPUNCH* IDATE*HD(60> 
COMMON  /LISTH/  OELTD  *  ALF AR , THRTL » G AMR  *  H I ND  *  VTF . WT , HF , DT , T I MS . XF  *FU 
lEL*MACH*VWF,LIFT*DRAG,THRST,RCF»ENGNO*THV« IN 
COMMON  /L I ST3/  CL  *  CD , S . OS ,TH ; R , ALTHR , CLMAX , AR  «  CL AR . CMACG . ALPHD ♦ CX  « 
1C2.0. IMOM*WINCR*CMDMP,OSC,ALMXR,VSF.BW 
COMMON  /LIST4/  FFS * JPOW « T GROS * DWA , THMOM • THREQ 
COMMON  /LIST5/  SIG,SOUND*NU*TEMR,PAMB* IATM,TEMF,DSODH.DRHO 
COMMON  /LIST6/  RTOD  *OTOR « KTOF  *FTOK ,NMTOF , FTONM  *RHOZ  *RH0Z2  *G2  *PZ  *  T2 
COMMON  /LIST8/  I NDEX  *C0EF  »GN»GT  *AN*AT  «  AX  *  AZ  * I TRM , VUPPER 
COMMON/L I ST9/ANS  *  ANS2 ,ND , X  *Y  ,Z  , 1 ND I C ♦ 1FLAS .NDER ( 3 ) « DBDX , DBDY , 

*  DBDZ«DB2DX,D82DY,0B2DZ 

COMMON  /LISTl 5/  IALMX, I NT , XEDOT . ZEDOT , XE . ZE , ! LOOP *DT I  ME , I  AT ,FL I FT  * 
1 DFLP*  C4HT .GALMXR 

COMMON  /LGEOM/  SHSW.SCR,OOOT ,RTH,DTHTH.RWA*DTHWA,DSTHZ,DWLH2, I Y Y , X 
1 CG . ZCG . CBAR «  PFN , S I  NTH , COSTH . ANGLE «  S 1 NAL  *  COS AL « RHZ ♦ DTHHZ . VTSQ .UDOT  * 
2W00T  *  THETR  *U«W 

COMMON/CONTROL/ JF I G . I RE V  *  I SP . NENG 
COMMON  /LIST99/  I ERR 

DEFINITION  OF  VARIABLES  IN  NAMELIST  • AERT 1  * 

S  *  WING  REFERENCE  AREA, 

BW  •  WING  SPAN 

C4HT«  HEIGHT  OF  THE  QUARTER  CHORD  ABOVE  THE  GROUND  (FT) 

GALMX*  ANGLE  OF  ATTACK  FOR  GROUND  CONTACT  (DEG) 

ALPMX  «  MAXIMUM  ALLOWABLE  ANGLE  OF  ATTACK  -  USUALLY  A  CL  LIMIT 
(DEG) 

ACLMD  ■  STALL  ANGLE  OF  ATTACK  (DEG) 

FPP  a  EQUIVALENT  FLAT  PLATE  AREA  OF  ADDITIONAL  DRAG  ITEMS  (SO  FT) 
WJNCO  ■  WING  INCIDENCE (DEG) 

DCLSP  «  INCREMENTAL  LIFT  COEFFICIENT  DUE  TO  LIFT  DUMPERS 
DCDSP  ■  INCREMENTAL  DRAG  COEFFICIENT  DUE  TO  LIFT  DUMPERS 
FODCL  *  INCREMENTAL  LIFT  COEFFICIENT  DUE  TO  ENG  OUT  CONTROLS 
EODCD  •  INCREMENTAL  DRAG  COEFFICIENT  DUE  TO  ENG  OUT  CONTROLS 
NX  *  NY  *  NZ  ARE  THE  NUMBER  OF  XC,YC*ZC 
XC  ■  THE  LIST  OF  ANGLES  OF  ATTACK  (DEG) 

YC  ■  THE  LIST  OF  THRUST  COEFFICIENTS 
ZC  •  THE  LIST  OF  FLAP  DEFLECTIONS  (DEG) 

CLA  *  LIFT  COEFFICIENT  AS  A  FUNCTION  OF  ( ALPHA « CT . FLAP  DEFLECTION) 
CPA  *  DRAG  COEFFICIENT  AS  A  FUNCTION  OE  ( ALPHA ♦ CT , FLAP  DEFLECTION) 
JF I G  *  1  FOR  MECHANICAL  FLAPS  PLUS  VECTORED  THRUST 
2  FOR  EXTERNALLY  BLOWN  FLAPS 


festers 
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C  3  FOR  INTERNALLY  BLOWN  FLAPS 

NAMFLIST  /AERT1/  S.BW. JFIG.C4HT .GALMX.  ALPMX, ACLMD.FBP. Wt NCD . 

1  NX . XC .NY . YC .NZ . ZC .CLA . COA , CCLSP.DCDSP .EODCL.EOOCO 

C  PFAD  INPUTS. 

tF  <  IN.NE.U  GO  TO  20 
READ  (CARD.AERT 1 ) 

WRITE  (PRINT. AERT1 ) 

WINCR»WINCD#DTOR 
GALMXR«GALMX*OTOR 
ALMXR*DTOR#ALPMX 
ALPHD  «  ACLMD»OTOR 
NO  •  1 
II  «  0 

COFP  »  FBP/S 
RFTURN 

20  RH0S*RH0Z2*SIG*S 

IF  (VTF  .LE .1.0)11*0 

OS  *  RHOS«VTF»VTF 
X  «  ACLMD 

IF  { I  I «EQ. 0  >  GO  TO  21 
IF( JFIG.LE.2)  GO  TO  23 
Y«TMMOM/QS 
GO  TO  22 
23  Y  *  T'.ROS/OS 
GO  TO  22 

21  Y  «  0 

22  7  *  OFLP 

CALL  LOOK (NX. XC .NY *  YC  «NZ  *  ZC .CL A .CD A ) 

CLMAX  «  ANS 

IF ( NENG.EQ.2 )  CLMAX«CLMAX+EODCL 
IF  (II.GT.O)  RETURN 
II  «  1 

VSF«SQRT(WT/(CLMAX#RHOS> ) 

RFTURN 

ENTRY  AER02 

60  ALTHR  ■  ALFAR-WINCR+THIR 
ALFAO  «  ALFAR4RT00 

IF(VTF.LE.l .O.OR. IREV.GT.O)  GO  TO  61 
tF( JFIG.LF.2)  GO  TO  65 
CT  n  THMOM/QS 
GO  TO  62 
65  CT*TGROS/GS 
GO  TO  62 

61  CT  »  O.OC 

62  X  «  ALFAD 
Y  «  CT 

7  *  OFLP 

CALL  LOOK ( NX . XC . N Y « YC . NZ . ZC . CL A . CD A  > 

IF ( JF IG.E0.2. AND . IREV.EQ.O )  GO  TO  72 

|F (VTF .LE. t . . AND. JFI G.FQ.3. AND. IREV.EQ. 0 )  GO  TO 

CL  ■  ANS 

cr>  »  ANS2 

I F ( ISP.GT.C)  GO  TO  66 
I F ( NENG.NE .2 )  GO  TO  64 
CL  ■  CL+EODCL 
CD  •  CD+EODCD 
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GO  TO  64 

66  CL  •  CL+DCLSP 
CO  *  CO+DCDSP 
64  CO  «  CO+COBP 
LIFT«CL*OS 
OPAG«CO*QS 
RETURN 

63  LIFT  *  THMOM+S I N  <  ALFAR+OFLP+DTOR ) 
OPAG  *  THMOM*COS<ALFAR+OFLP*DTOR) 

rfturn 

72  CL  »  ANS-CT#SIN(ALTHR) 

IF  < VTF  .LE.l.O)  GO  TO  73 

CO  ■  ANS2+CT#C0S(ALTHR) 

GO  TO  74 

73  CO  «  ANS2+CT*C0S<ALTHR) 

74  CO  «  CD+CDSP 

IF  <|SP.GT*C)  GO  TO  75 
IF  (NENG.NE.2)  GO  TO  80 
CL  »  CL+EODCL 
CD  ■  CO+EODCD 


GO  TO  80 

75  CL  ■  CL+DCLSP 
CO  *  CD+OCDSP 
80  LIFT •CL+QS 
rRAG«CO*OS 
RFVURN 
FNP 


P. 


FUNCTIONAL IN( I *X* Y*XA) 

DIMENSION  XlthYtl) 

alin  ■  y<  i  >+(xa-x(  m#(v<  i+n-vm  >/<xu+n-x<  1 1 > 

RETURN 

FNR 
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SUBROUTINE  AT MO 5 

C  6400  FORTRAN  ATMOSPHERIC  SUBROUTINE  WITH  MJL  STO  210A  TEMPERATURE  OPTIONS* 

c  this  subroutine  calls  for  the  following  entries. 

C  KABO « OK ABO 

C  INDEX  IATM«  SET  IN  THE  MAIN  PROGRAM.  SELECTS  THE  TEMPERATURE  OPTION. 

C  I  ATM  «  0,  FOR  STANDARD  TEMPERATURES 

C  I  ATM  *  I  MIL  STO  2lCA  TROPtC  TEMPERATURES. 

C  I  ATM  a  2  MIL  STD  210A  POLAR  TEMPERATURES. 

C  I  ATM  a  3  MIL  STD  210A  HOT  TEMPERATURES. 

C  I  ATM  a  4  MIL  STD  210A  COLD  TEMPERATURES. 

C  I ATM  a  5  TEMPERATURE  IS  INPUT  AS  TEMF 

INTEGER  CARD . PR ! NT » PAGE 
PFAL  LIFT. NU. MACH 

DIMENSION  TROPICUS)  .P0LARI20)  ,H0T(20)  .C0LD<35> 

COMMON  A.IST1/L INE.PAGE.LIMIT.CARD.PRINT. INP. 1 PUNCH, I DATE. HOC 60 > 

COMMON  /L I ST2/0ELT0 • ALFAR , THRTL » G AMR , H  J  NO « VTF , WT . HF , OT • T 1 MS . XF , 

#  FUEL . MACH . VWF • L I FT .DRAG , THRST . RCF , ENGNO « THV , l N 

COMMON  /L I  STS/S  I G . SOUND . NU . TEMP , PAMB , I  ATM « TEMF . DSODH  «  DRMC 
DATA  NTROPC.NPOLAR.NHOT.NCOLD/?  .4  «4  *  7/ *  ( TROPtC  (  1  >  .1*1 ,15)/ 

1  -84270. 12.62808.804,546. 14631 .0039108943,53595. .-2J5BC3.01. 

2  94397,625.223.43864, «0C2 03 3 033^.69620. .453.66956,67598.227. 

3  289.9406.0,0013433339. lOCCOO./. (POLAR! I ) . I  a J , 20 ) /220 . 31 1 53 • 

4  -1058,4386,  443.77985,  C.00310U192.  3243,,  1623.9492.  15111.866. 

5  455.6965.  -5 . 1 608824E-C4 .  9882..  15029.377.  39599.146.  478.47166. 

6  -0,002800841,  30065..  3663.7771,  13319.676,  400.08921. 

7  -2.5384026E-O4,  86C92./. (HOTU ) » 1 =1 .20)/- 196560. 19.  55620.773. 

8  559.154C7,  -0.0039742099,  39400 1 09942 ,52 .  73731,906. 

9  390.98426,  3 . 1 72946E-04 .  504C0. .- 1 23 .23365 ,  47876.697.  410.41972. 

#  I  .8486342E-0A,  66400.,  10380998.,  -63820.346.  260.37596. 

1  O. 0012429502,  i',0C00./.(C0L0(!  ) « I  a  1  »  35  > /-4 8393 »  349,  10737.404. 

2  395,18101.  0.012984173,  3311.,  0,0.  1.0.  444.688.  0.0.  1C744.. 

3  627229540.,  381783.16,  2121.0792,  0.00131C1285.  30715*.  0.0.  1.0. 

4  374.688,  0.0,  42377,, -1609. 4722,  51071,014,  597,1)641, 

5  -0.0052531685,  50583.,  0.0,  l.C,  334.688.  0.0,  61007., 

6  -5071005 0,.  -31326.366,  1046,3219.  -0.0026667396.  73055./ 

N AMEL l ST / A TM 1 /NTWQPC  • TROP 1 C . NPOL  AR . POL  AR , NhOT  «  HOT  « NCOLD , C OLD 

1F< IN-1 )5» 1 .5 
1  WPf TE(P9!NT,ATVJ ) 

Oc-TURN 

C  STANOARO  atmosphere 

5  IF(MF. OF. 36089,239)00  TO  10 

PAMB  «  21  1 6,221  8*  (  I  .-,687586  “>64?-05*kF  )  a*5.2S58761 

IFUATm,GT.O)GO  TO  15 

Tr«R  *  1 .8*(288,I5-.CC19«12*HF) 

GO  TO  30 

10  IFJHF.GE. 65616. 798)00  TO  11 

PA**P  »  472,60)64»EXPU.7J4S725-,480634?PE-;54*mE) 

IF( 1 ATm.OT ,0 ) GO  TO  15 


fs  . 

I 


y- 

E; 

i; 


ft 

r 


I 
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Trva  •  359.97 

60  TO  30 

t)  PAMfl  •  1 14.3457/1. 9076852+.14068?74E-OS*HF)**34.  (63194 
IFltATM.6T.0lG0  TO  IS 
TFMR  >  I *6#t I96.65+»0003048*HF) 

60  TO  30 

15  60  TOt 16*1000.2000.95.261 « 1ATM 

C  TROPICAL  ATMOSPHERE  TEMPERATURES. 

16  IFtHF.GT. 100000. 160  TO  25 

CALL  KA0O<NTROPC*HF«TEMO. TROPIC! 

60  TO  30 

25  TFMO  »  424.27 
60  TO  30 

C  POLAR  ATMOSPHERE  TEMPERATURES. 

iOOO  lFtHF.GE.86092.!GO  TO  1130 

CALL  KABOtNPOLAR.HF. TEMP. POLAR) 

60  TO  30 

1130  TFMR  ■  378.27 
60  TO  30 

C  MOT  ATMOSPHERE  TEMPERATURES. 

2000  IFtHF.GT. 100000. )60  TO  2120 
CALL  KA0OtWOT.MF.TEMR. HOT) 

GO  TO  30 

2120  TEMR  *  448.07 
60  TO  30 

C  COLO  ATMOSPHERE  TEMPERATURES. 

95  IFtHF.GT. 100000.  )60  TO  170 
CALL  KA0DtNCOLO.Mr.TEMe.COLO) 

GO  TO  30 

170  TEMP  a  355.77 
GO  TO  30 

26  Tr*B  a  TEme.«39,*7 

30  516  •  ,245C0246*PAM8/TEMR 
TF-F  •  TEMR-459.67 
TfMK.TEKRyi ,Q 
SOUNO»65.77a35C*SOOT  t  TE  -wC  | 

MJ  *  .1261 ‘25E-044SORT <  TEMK }/<SIG*< 1  ,♦! 1Q.4/TEHKI I 

of turn 

mho 


a 


1 


1 

I 


V 
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SUBROUTINE  P  I  NO  ( N  i  X  •  Y  »XC  *  YC ) 

COC  6400  SUBROUTINE  WHICH  PERFORMS  A  LAGRANI AN  INTERPOLATION 
ON  A  ONE  DIMENSIONAL  ARRAY,  IF  THE  ARGUMENT  IS  OUTSIOE  OF  THE 

RANGE*  A  linear  extrapolation  is  maoe* 

this  SUBROUTINE  CALLS  the  FOLLOWING  ENTRIES 
INTP2 

N-  THE  NUMBER  OF  X»Y  POINTS  IN  THE  ARRAY 

X-  THE  ARGUMENT  FOR  INTERPOLATING 

Y-  THE  SOLUTION 

XC-  THE  ARGUMENT  LIST 

YC-  THE  DEPENDENT  VARIABLE  LIST 

DIMENSION  XC(20) «YC(?C| 

Integer  card* print, page 

COMMON/L I ST I /L I NE *P AGE *L 1 M 1 T • CARD *PR I NT  * i NP« I PLMCH , I DATE • HD ( 60 ) 

1FfN.LT.4l  GO  TO  SC 
l»l 

N|*N-I 

00  10  I*| »NI 

1FfX.GE.XCt I ! * AND.X.LC .XC ( 1*1 II  GO  TO  IS 
|0  CONTINUE 
GO  TO  90 

19  IFtX.NB.XEf | »»  GO  TO  20 

V*vc«|! 

RETURN 

20  IFtx.NE.XCt  l*n»  GO  TO  25 

v*vctm» 

RfTURN 

25  |P( I.GT.2I  GO  TO  30 
1*1 

GO  TO  49 

30  tFtf.LY.N~l!  GO  TO  35 

l*N~3 

GO  TO  45 
35  t*f-l 

49  NX* 1*3 

CALL  !NTP2lNX.Xe«YC*S«V,A,R»C,0l 
of turn 

50  IFtX.LY.XCI  I  I >  1*1 
1FtX.GT.XCtN!!  t*NI 
V  *  ALINtt,XC*YC*X| 

RETURN 

ENO 
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SUBROUTING  GILL( tuOT «V«OUMO« ILOOP tDTJMEj 

THIS  IS  AN  INTEGRATION  SUBROUTING.  TO  EXECUTE  AN  INTEGRTATICN  STE 
OEOUIRES  FOUR  PASSES  WITH  I LOOP  SET  TO  2  THROUGH  5.  THE  FIRST  TIME 

subroutine  is  used*  the  first  pass  is  made  with  iloop  set  eoual  to 

ALL  VARIABLES  BEING  INTEGRATED  MUST  BE  CALLED  IN  EACH  PASS  BEFORE 
GOING  ON  TO  THE  NEXT  PASS. 


PV«VOOT*OT!ME 

o*nuMO 

GO  TO  I l 0.20 .30 *40 #50) * 

10  0*0. 

C*l ./SQRTI2. f 
CMINUSatfC 
CPLUS*1.+C 
20  e».5*0V-0 

0*3,«n-.S«0V+Q 
GO  TO  60 

30  »«pP!NUS*«DV~0| 

0*0*3 • «R-CM INUS»0V 
GO  TO  60 

40  O*CPLUS«f0V-0j 
0*<H3.«R~CPLUS«0V 
60  TO  60 

so  e*<oy*2.*0)/6. 

0«0+3*»e-.5«OY 
60  V*V*R 
OUMO*C 
OfTuBN 

eno 


iloop 


u 

u 

u 

u 

u 

if 

u 

u 

u 

u 

u 

u 

u 

u 

u 

u 

u 

u 

u 

u 

u 

u 

w 

u 

u 

u 

u 

it 


■m 


m 


■fi 
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SUBROUTINE  6INT6  A| 

C  6400  FORTRAN  SUBROUTINE  WHICH  INTEGRATES  At.TITUOE.R/C  AND  ATTITUOE  At 
C  AND/OR  6ROUNO  SPEED*  DISTANCE  At 

C  OR  GROUND  SPEED.  01  STANCE  AND  ATT1TU0E  At 

At 

c  this  program  calls  the  following  subroutines*  ai 

C  AEQFM, ATMOS, GILL  A| 

At 

INTEGER  CARO. PO| NT* PAGE  At 

REAL  LIFT.KTOF.MACH.Ng.NMTOF.IYV  At 

COMMON  /LIST!/  LINE. PAGE. LIMIT. CARD. PRINT. INP. IPUNCH. | DAT E . KO ( fcO J  AI 
COMMON  /LISTS/  OELTD. ALFAR.THHTL .GAMR.hlNO. VTF  »  WT .MF.OT.TIMS.XF «FU  AI 
JFL.MACK.VWF .LIFT. DRAG. ThRST.RCF.ENGNO.THV. IN  A| 

COMMON  /L 1ST J/  cl.co.s.os.thjr.altmr.clmax.ar.clar.cmacg.alpho.cx.  AI 

It  .  |N;_».C  D  iP.OSC.ALMXR.VSF.BW  A| 

COMMON  /lISTA/  FFS.JPOW.TGROS.OWA.THMOM.THREO  A I 

COMMON  /LISTS/  SIG. SOUND. NU.TEMR. PAMS, |ATM. TEMP, OSOOH.QRhO  Af 

COMMON  /L 1 ST6/  PTOO.DTOR.KTOF.FTOK.NMTOF.FTONM.RH02.RH02S.G2.P2.r2  AI 

COMMON  /LISTS/  |NCEX, COfF.GN.GT.AN, AT, AXE.A2E. ITRM.VUPPER  A| 


COMMON  /LIST  1 5/  I AlMX* |KT .XEDOT .2ED0T .XE.2E . ILOCP.OTIME, I  AT .FlIFT.  At 
t OFLR . C4HT • G ALMXP 

COMMON  /lGEOM/  SHSW. SCR. BOOT. RTH«0TMTH.RWA.DTHWA.0STH2.0WLH2*tYY,X  A| 
t CG.ZCG.CSAR.PFN. StNTH. COSTH. ANGLE. StNAL.COSAL.RH2.OTHM2.VTSO.UDOT .  A| 


SWpOT.THETR.U.W  AI 

AI 

10  IF  IXEOOTt  30.20.30  A I 

26  6AM0»O«O  At 

GO  TO  40  A| 

30  GAMR.ATANUEOOT/XEOOTj  A| 

AO  VTS0*XED0T«XF0oT^2eDOT*2eO0T  A I 

VTF.SQRTtVTSOl  AI 

alfar«thetr-gamo*wincr  At 

IF  IfALMXt  60. 60. SO  At 

90  aufar-almxr  A| 

0« 1 xEOOT*2E20T-2eOOT#XE2UT  t/VTSG  A  t 

60  CALL  ATMOS  At 

MACM»VTF/50UN0  At 

CtS*RM022*SIG4S«VTS0  AI 

CALL  AEOFM  A I 

XF20T-AXE  At 

»*-20T*a2f  AI 

CO  to  f TO. TO. SO, TO, @01,  IlOQP  AS 

TO  T|Mi|«ttMS*.9*0TtMF  Af 

SO  XC/iOT •  X?D$T»VWFF  AS 

CALL  GILMXGOOT. XF.OSff.  SLOOP. OTtwet  AI 

CALL  6tLLtXE20T.XEtOT.OxE2.IL<H)P.OTI^;t  At 

t *  eiNT~2»  IIO.IOO.RG  AI 

DO  Call  GILLt?«OOT,2E,OCe.  ILOOP.OTImsi-  I  At 

CALL  StLLtge20T.2€oOt.O*£2»tLDDP.CTJMe>  At 

too  Call  e»LL«a.?*KfT#.arH.iLOoe.GT|‘*£»  as 

tto  GO  TO  £125.530. |>C.t3r.tACj.  SLOOP  *1 

120  ILOOPA3  A I 

60  TO  10  *! 

ISO  tLDOP*ILOOP#S  A I 

GO  T0  to  AJ 

540  ILO-0P*2  A I 

pf-TuRN  .At 

End  At 
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SUBROUTINE  HEAO(NL)  AK 

AK 

INTEGER  CARQ.PR JNT .PAGE  AK 

DIMENSION  IDI  I C  |  »KP(  IC)  »N{_  <8 )  AK 

COMMON  A.ISTI/  LIN€»PAGE.UHIT,CAR0*PR1NT*|NP»IPUNCM.  J0AT£«H0C60>  AK 

AK 

00  2**  I»I*4  AK 

AK 

IF  <MOO(I.2).N£,0)  WRITE  {PRINT. 30)  AK 

30  FOR* AT  (1H  )  AK 

PFCO0£  < 10.«0.NCt Ji )  10  AK 

*0  FORMAT  t«l 2)  AK 

pO  JO  K« 1  *  10  AK 

KK«JOtK)  AK 

10  MP<Kt*KDC<K)  AK 

WRITE  (PRINT. SO)  HP  AK 

SO  FORMAT  < I2X.10IEK.A10) )  AK 

20  CONTINUE  AK 

LInR*UINE+6  AK 

OFTURN  AK 

FNft  AK 


SUBROUTINE  l-!TPJMM>«»X.V«Xc*YC*A.e.C.O> 

C  6400  FORTRAN  3Rf)  0F6»EF  POuVNOMUt  1 NTERPCU A T l ON  SUBROUTINE, 

C  PNTOV  with  JluasCRiPTFO  x  A  NO  V  Q6TA. 

DIMENSION  XI) J«  V(|| 

S  K  NK-? 

10  <*U  *  {¥ t J )-Y( $«) ) »/(*( 1 |~XI 1*1 > | 

&C  *  (AU-tTI  })•*¥(  U2l)/|K>  |  )-Xt  !*2)  I  »X»Xt  1*}  J-X(  !*gt  } 

A  s  tfff -««(.-(¥(  1  |-Vt|*3»)/0»<!  }“X(  lOHI/iXIU)  )-X(  iOU)/iXtt*£)o 
•X( t*3)  l 

na  »  x  1 1 ) *s  f | » t  } 

P  *  *t£«A.#(^A*X(  1«?|  I 

0  a  »1,«A4  txi  t  1  *X(  I  l«M(  1  •  1  )*^A)«{»*^* 

0  «  T(  f  J*f  ti(Si|  |  j  .‘•(♦S'l  |  )*€)*»!  I  ) 

VC  *  » (*«»XC**)*xt><l*XC*0 

Cfc-'TyRN 

F*«0 
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SUBROUTINE  KABDINF1TS.XC.VC.X| 

C  6400  FORTRAN  KABO  EVALUATION  PROGRAM, 

C  NT  ITS  IS  THE  NUMBER  OF  CURVE  FITS  STRUNG  TOGETHER  FOR  THE 

C  PARTICULAR  VARIABLE. 

INTEGER  CARO.PRlNT.PAGE 

COMMON  /LIST!/  LINE .PACE .LIMIT .CARO. PR I NT. INP. I PUNCH. 1 DATE. HO (60) 
COMMON  /LIST99/  I ERR 

DIMENSION  X(50)  X 

IfRRvO  X 

1*1  X 

10  IF  (XC«X<|A4tt  ?0 . 20 . 30  X 

20  Ve»XII»X<XC-X||AI»|AKIlA2|AXIIO>*XC  X 

OFTUBN  X 

SO  1*1*5  X 

IF  (I.LE.SANFITS!  GO  TO  10  X 

WRITE  (PR1NT.40)  NFITS.XC.Xt  t-U  X 

IFRR*  I  X 

OPTURN  X 

C  X 

40  FORMAT  (SX.3SM  UPPER  LIMIT  OF  X  EXCEEDED.  NFITS*. I3.3H.  XC-.E12.3  X 
I.9H.  XI !♦«)•. E12.91  X 

f'Nft  X 


AuPROuriMI  LAGOfllM.Y. |NO*x>  Y 

C  64CU  FORTRAN  StUBSGU  T I  Ng  *N|CM  OOCS  A  FOUR  POINT  VAGRANGIAN  INTLRPO  V 

V 

OIMEHtllON  XtSl.ytOt  .OXt  |U|  v 

V 

tr  (tNCZvy  33.|0. SO  V 

»0  1  1*1  NOCK*  S  Y 

Y 

t>0  1*1.4  Y 

*  I*  t*  1  V 

tv*  20  3*11 .5  Y 

r*X*I  Y 

Y 

to  CQNTjNruC  Y 
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SUBROUTINE  LANDING 


C  6400  FORTRAN  IV  LANDING  SUBROUTINE  FOR  5T0L  AIRCRAFT 

C  TWF  BOUNDARIES  OF  THIS  PROGRAM  ABF  - 

C  VTD  GT  OR  £0  TO  APV«VSTALL 

C  V08S  GT  OR  EO  TO  APV*VSTALL 

C  TOUCHDOWN  RATE  OF  S|N<  LT  OR  EO  TO  RS 

C  THETA  AT  TOUCHDOWN  GT  OR  £0  TO  ZERO 

C  THETA  AT  TOUCHDOWN  LT  OR  EO  TO  GALMX  |A£RO  INPUT) 

C  ALPHA  DURING  APPw.dCM  LT  OR  EO  TO  ALPMX  (AERO  INPUT) 

c  this  program  call  thc  following  subroutines 

C  ATMOS • SX I P»H£ AD. TMST 1 .NWTP2.AER01 .GINTG 

INTEGER  CARO. PRINT, PAGE 
OFAL  LIFT. <TOF. MACH, NU.NMTOF 

OIMCNSION  NlO) 

COMMON  /LIST!/  LINE.P^GE. LIMIT, CARO. PRINT, INP. IPuNCH. I0ATE.HOI60I  F 
COMMON  /CIST2/  0£L7C.AlFAR,TWHTL.GAMR,mINO.VTF,wT.HF«OT,T|MS.XF,Fu 
IfL  *  MACH . VWP .LIFT. ORAG  «  ThRST , RCF . EN&NO , THV . I N 
COMMON  /LISTS/  CL«CD,S.a5,YHtR,ALTMR,CLMAX.AR.CLA».C>«ACG.ALPMD»CX» 
IC2.0. IMOM.W'NCR.CMOMP.OSC.ALMXR.VSF.aw 
COMMON  /LISTA/  FFS. JP0W.T6RGS.0W a  9 TKMOM.TmREO 
COMMON  A. t  STS/  «|G, SOUND. NU.TEMR.PaW9, | ATM . TEMF .QSOOH.ORhO 
COMMON  /LIST6/  RTOO.DTOH.<TOF,FTOK.NVTOF.FTONM,RmOZ,RmOZ2,GZ.PZ.TZ 
COMMON  /LISTS/  lNOex.COFF.GN.  jiT.AN.AT.AX.A2,  ITR.V.VUPPER 
COMMON  /LlStiS/  lALKX.lNT.XEOOT.ZEOOT.XE.Ze.lLOOP.OTlMe.lAT.FLIFT. 
|CFlO.C*hT.GALMXB 

COMMON  /lGFOM/  SHSW.SCR.OOOT .RTH.CTHTH.RWA.orMWA.DSTHZ.OWLHZ, IVV.X 
l CG. 2CG.C9AR .BFN, SI nTm, COSTh. ANGLE , S I NAL , COS  AL  »  Rh£ .OTmmZ . VTSO.UOOT . 
JWfVJT  *  THE TR ,u»W 

COMV0N/C  ONTROL / jf t  L . I RE  V . I SP . NCNG 
CO***»on  /LIST99/  IEOO 

PA*  A  nFLP.R0TPT.HF0ft.R0TATN.T0OX  « YSP , T»€ V , APR . SCOEF »aCO€F .DGTm . 
)rnTn./6C'».i,.Fo.,0.,?.,2,,?.,  t ,)  ,o,: .o.a.c.T.tuis/ 
oat *  gamm* ,os/*t.m( io*/ 

OATa  >  NL  m  .  |  •  I  ,3  )/ 1  nN'“S*j6TT02S6.  S  .  ISK393SJSI  *<5*2  .  i 

1*OOISISI  *  I0»0»6tct  l0H3Ot09t  |$8,  )  CMS?3R>V3e:><3.  »  CmS9*0&»5I  30/ 

€  DEFINITION  QF  vaRIAHlES  in  nam£lIST  *l*noi* 

C  OFlR-  Flap  SfTT|N4  tOEGi 

c  obstacle  height  »ft> 

C  O0TATN*  ROTATION  RATE  (DEG /SEC  > 

c  v«£sx*  minimum  control  ..««  eg  on  the  attend}  j«ts> 

C  VMCAK*  minimum  CONTROL  SPIED  IN  THE  AIR  t  *ST  S ) 

C  YfHJK*  T|MF  OfL  Ay  AFTCP  roue mSDaN  FOR  OLAKE  APPLICATION  t  SEC  I 

t  Tjp»  T|xf  Of  LAV  AFTER  TOUCHDOWN  r©«  JPOlL*®  SCPlOYMEnT  4  SEC  < 

r.  TRFVA  TIME  OfL Av  *FT«»  TCuC«00*N  '0»  THREAT  RfvtBSAL  tSFCl 

e  ARC*  Ratio  of  aP’PcoaCh  S^CEO  to  minimum  cgiNTROL  SPEED 

C  OGTH"  | NCRF MfNit *L  LOAD  FACTOR  AVAILABLE  AT  Th?  OBSTACLE 

C  v?GtO«  |NC«^MfW»AL  LC*y  FAvTCil  AVAILABLE  AT  TOUCHDOWN 

t  CPCOfF*  SKILL  *  NO  CCEFFICICNT  OF  FRICTION 


C  BCOEF*  BRAKING  COEFFICIENT  OF  FRICTION 

C  GAMMA*  INITIAL  FLIGHT  PATH  ANGLE  IOEG)  (NEGATIVE  IS  DESCENDING) 

C  OS*  MAXIMUM  RATE  OF  SINK  AT  TOUCHDOWN  (FT/SEC) 

C  (POSITIVE  IS  DESCEND  NG 

C  OOTPT*  EQUAL  TO  ZERO  SUPRESSES  PRINTING  OF  THE  ENTIRE  TIME  HISTORY 

C  ThF  FOLLOWING  VARIABLES  ARE  ENTERED  at  TIME  OF  LOADING  AND  ARE  USED 

C  UNTIL  OVERBIDDEN  BY  REAOINC  THE  APPROPRIATE  VARIABLES  IN  LAND! 

C  r*LP  ■  60.0 

C  HFOtt  ■  50.0 

C  ROTATN  a  B.O 

C  T8RK  ■  2.0 

C  TSP  »  2.0 

C  TREV  •  2.0 

C  APR  •  I .10 

C  OGTH  ■  0.30 

C  OGTO  •  O.IS 

C  RCOEF  *0.10 

C  BCOEF  •  0.30 

C  gamma  ■  -T.90 

e  rs  •  io. o 

C  ROTPT  »  1.0 

NAM£L 1 5T/L AND  I /OFLP . »0 tPT . MF 08 . ROT ATN . VMCGK » T0RK . TSP • APR . RCOEF ,B " 0 

1 FF. TREV .GAMMA .RS . VMCAK . OGTH. OGTO 

IF  ( IN.NF.I)  GO  TO  5 
Pf AO( CARO. LAND) ) 

WRITE (PR I NT .LAND) ) 

OF TURN 

5  mfbun.mf 

mf«hFRun>c*hT*3». 

CALL  ATMOS 
OmOS*RmOZ2*SIG*S 
V«CG  •  VMCGK**CTQF 
VMC  A * VMC AK*K  T  OF 
APVaAPR 
CO*F  •  0.00 
|-5P*0 
JCa-VsO 
•  O 

JWV  *  0 
j*e*  •  0 

•  o 

WRITE  («H3|nT,  I 

(CCO  FORMATl  )  M)  .3;«,  *(5uTOyT  DEFINITIONS  -  LAND)  NO*  .  /  •  «R  «  *5EGM€NT  •  *, 

1*1  CETERMINAT  ION  of  ancle  of  ATTACK  AT  TQUCMUQWK* • / .  |*J£. 

?•?  R/C  AVAtLAPu?  AT  TOuCkDOwN  -  \  ENGINE  OuT*./,J*K, 

3*S  DETERMINATION  OF  ancle  OF  ATTACK  AT  the  OC’STAClE*./. I*x, 

«•«  R/c  AVAILABLE  AT  Twf  OBSTACLE  -  5  ENGINE  OuT*./,l*X, 

9*B  INTEGRATION  P80N  OBSTACLE  TO  TOUChDOWn*. /» >**, 

/**  YouCmOOWN  to  NO?FCKJW«  INTEGRATION*./, I «K, 

?**?  CONDITION*  at  NOSE  DOWN*./.  I  AW. 

«*«*  NOSFOO'WV  TO  stop  INTEGRATION*. /» |«X, 

R*R  CONDITIONS  at  STOP*,/) 
line  •  L«K»T*l 
GO  5C 

C  OUTPUT  BLOCK 


20  IF  (R0TPT.NE.0, )  GO  TO  2.5 

IF( (ISFG.E0,6).0R.( ISEG.EQ.8) )  GO  TO  35 
25  IF  (LINE.LT. LIMIT)  GO  TO  30 
CALL  SKIP 
CALL  HEAO(NL) 

WRITE  (PRINT. IOOII  H0(27) 

1001  FORMAT  (1H+.IX.A10) 

30  VTK*VTF*FTOK 

GAMD»GAMB*RTOD 
OCM«60 ,*VTF*S I N ( GAMR ) 

ThETD»THETR*RTOD 

thoto»q#rtoo 

ALFAD»A!.FAR#RT0D 

althd«althr#rtoo 

WRITE  ( PR  I NT . 1 002 )  I SEG  v HF . DFLP . WT , VTK . RCM , THETD . AN .CL .L I FT . TGROS . 
1TIMS.XE.ZE.GAMD*  ALFAD*  THOTD. AT .CO . DRAG. ALTHO 
L  t NF«L  t  NE+3 

1002  FORMAT  ( 1H0. I  11 .6F12.3.2F12.5.2F12.3./. 1 2X .6F1 2.3 .2F12.5.2F 12.3 ) 

(THRTL.GT. I .0)  GO  TO  AC 
IF  (THRTL.LT .0,0 )  GO  TO  45 

35  GO  TO  (125.150.1 75*200 .250.260.300.330. 350 ) .  ISEG 
40  1F( ISEG. NE.l .OR. ISEGoNE. 3. OR. ISEG.NE.5)  GO  TO  35 
WRITE  (PRINT. 1008) 

1008  FORMAT  (1H  «13X,#THR0TTLE  GREATER  THAN  MAXIMUM  POWER*) 

GO  TO  35 

45  IF( ISEC.NE.l .CR. ISEG. NE. 3. OR. ISEG.NE.5)  GO  TO  35 
WRITE (PRINT .IGOR) 

1009  rORMAT ( 1H  . 1 3X ,*THROTYLE  LESS  THAN  IDLE  POWER*) 

GO  TO  35 

C  SFY  APPROACH  SPEED  -  FACTOR  TIMES  STALL  SPEED  (POWER  ON  1  ENG  OUT) 

50  lN*0 

ThRTL«1 .0 
CALL  AERO  I 

ALTMR« ALPHO-W I NCR+TH IR 
VTF1«VSF 
ENGNO-ENGNO-1 . 

NFNG«2 

MACH- VTF 1 /SOUND 
VTF«'-'TFI 
NSTFP^O 
55  CALL  THSTl 
CALL  AER01 

IF ( JF I G.EQ.2. AND.VTF.GT • 1 • )  CLMAX«CLMAX- ( { TGR0S*S IN ( ALTMR ) )/ (RHOS* 
1VTF*VTF) > 

VTF«SQRT( (WT-TGR0S*SIN( ALTMR) ) /<RH0S*CLMAX ) ) 

YVTF« ( VTF1-VTF ) /VTF 
IF  (ABS(YVTF) ,LF. 0.0001 )  GO  TO  65 
IF  (NSTEP.EQ.O)  DMAC*  MACH-VTF/SOUND 
CALL  NWRP2( DMAC. YVTF, MACH, YSAV.XSAV) 

NSTEP«NSTEP+1 
IF  (NSTEP.GT .15)  GO  TO  60 
VTF 1 «MACH*SOUNO 
VTF»VTFl 
GO  TO  55 

00  VTF »MACH*SOUND 
65  VSF-VTF 

ENGNO-ENGNO+1 . 

NFN6»0 

!N«2 


VSX»VSF*FTOK 

VAPF«VMCA*APV 

IFCVAPF.LT.VMCG*APV)  VAPF«VMCG*APV 
IF(VAPF.LT.VSF«APV)  VAPF«VSF*APV 

IF{ VAPF.LT.VSF*SORT< 1 .+OGTH) )  VAPF«VSF*SQRT t  1 .+DGTH) 

I F ( VAPF ,LT , VSF#SOST  < 1 . +OGTD ) >  VAPF ■ VSF*SORT  C 1 • +DGTD  > 

APV  ■  VAPF/VSF 
t  APV*0 
IGAM  «  0 
70  VTF«VAPF 

V6PK«VAPF*FT0< 

DFTFRMJNF  ANGLE  OF  ATTACK  AT  TOUCHDOWN 

100  VTF«VAPF 

TIMS  «  0.00 
xe«o.o 

2F  «  0  »  0 
0»0,0 
InTEG*0 
NF«HFRUN+C4HT 

IF  (tGAM.EQ.O*  OR .GAMP «  LT  » G AMM A*OT OR )  GAMR»GAMMA*OTOR 

CALL  ATMOS 

RH0S«RH0Z2*S I G#S 

0S«RH0S*VTF#VTF 

MACH«VTF/SOUND 

1 TRM  «  1 

GN*1 .0 

GTkO.O 

lNOEX«6 

C0EF-0.0 

CALL  AEOFM 

THETR« ALFAR+GAMR- W I NCR 
IF  CALFAR.LE.ALMXR)  GO  TO  HO 

INCREASE  VELOCITY  TO  DECREASE  APPROACH  ANGLE  OF  ATTACK 

105  t  APV« I APV+1 

APVl«APV+FLOAT< IAPV1/100, 

VAPF'»VSF*APV1 

WRITE  (PRINT. 1003)  APV1 

1003  FORMAT < IH  .13X.1&HVAP  INCREASED  T0.F7.3. 12HTIMES  VSTALL*/) 

LINE  «  LINE+2 

GO  TO  70 

DECREASE  GAMR  TO  OECREASE  R/S  TO  MAXIMUM  INPUT  VALUE 

110  tF  <VTF*SlNCGAMR>.GE.-(RS*i.Oi ))  GO  TO  115 
GAMR«A5 INC-RS/VTF ) 

igam»igam+i 

IF  CIGAM.GT.20)  GO  TO  405 

GAMO  »  GAMP4RT0D 

WRITE  (PRINT. 1004)  GAMD 

1004  FORMAT ( 1H  . 13X.*MAXIMUM  RATE-OF-S I NK  EXCEEDED.  GAMMA  RESET  TO* 
1F7.3.*  DEG*) 

LINE  ■  LINE+2 
GO  TO  100 

CHFCK  THETA  TO  AVOID  ‘TAIL  HIT*  OR  ‘NOSEWHEEL  FIRST*  LANDING 


115  IF  (THETR.GT.GALMXR-WINCR)GO  TO  105 


o  o  o 


;wwssa«*s!so*K3s»»wsas3«a«BBSf«s‘*Bowi«<*s>nB!racia*»^^ 


'/' .  --‘2Ci^r®^swwo'aR»OiS3WW^a>a«»xiwcap.«a'<i5i 


IFlTHETR.GE.-O.OOi  >  GO  TO  120 

GAMRaW l NCR-ALF AR 

IGAMalGAM+l 

IF  U6AMcGT.20-,  GO  TO  405 

IE<GAMR«LT.0.<  ,AND.VTF*SIN(GAMR>.GE.-RS>  GO  TO  100 
GO  TO  400 
120  15FG  «  l 

GAMTR  a  GAMR 
60  TO  20 

C  DETERMINE  R/C  AVAILABLE  AT  TOUCHDOWN 

125  THRTL  »  1.0 
index  «  3 

ENGNO  a  ENGNO-l.O 

NENG»2 

CALL  AEOFM 

ENGNO  »  ENGNO +1 ,C 

NENG*0 

THFTR  a  ALFAR-W 1 NCR+GAMR 
ISE6  *  2 
GO  TO  20 

C  DETERMINE  ANGLE  OF  ATTACK  AT  THE  OBSTACLE 

150  HE*HFRUN+C4NT+HF08 
ZE  a  HFOB 
GAMR  a  GAMTR 
CALL  ATMOS 
RH0S«RH0Z2*$ I G*S 
OSaRHOS»VTFaVTE 
MACH-VTF/SOUND 
GN«1 .0 
GT»0.0 
INDEX  «  6 
ITRM  a  1 
COEFaO.O 
CALL  AEQFM 

THETR«  ALFAR+GAMR— W I  NCR 

IF  < ALFAR.GT.ALMXR)  GO  TO  105 

I SEG*3 

V08S-VTF 

GO  TO  20 

C  DETERMINE  R/C  AVAILABLE  AT  THE  OBSTACLE 

175  THRJL  ■  1.0 
tNDEX  a  3 
ENGNO  a  ENGNO-l.O 
NFNG-2 
CALL  AEQFM 
ENGNO  a  ENGN0+1.0 
NENGaO 

THETR  -  ALFAR-W I NCR+GAMR 
I ?EG  »  4 
GO  TO  20 

INTEGRATION  FROM  OBSTACLE  TO  TOUCHDOWN  (DUE  TO  THE  NO  FLARE 
LANDING  TECHNIQUE  USED  IN  THIS  PROGRAM  THE  INTEGRATION  IS  MADE 
IN  ONE  STEP) 
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200  GAMR  «  GAMTR 

IFIGAMR.GE.0.01  GO  TO  400 
IF  ITHETR.LT. -0,001 >  GO  TO  ACS 
TIMF1«HF0B/(VAPF#SIN(ABS<GAMR} ) ) 

P1ST1-  T I  ME  1 *VAPF#COS<  GAMR ) 

TIMS-TIME1 

XF-DIST1 

ZE-0.0 

HF  ■  HFRUN+C4HT 

CALL  ATMOS 

RH0S«RH0Z2*S I G#S 

OS»RHOS*VTF*VTF 

ITRM  «  1 

GN-t .0 

6T-0.0 

INDEX-6 

C0FF«0*0 

CALL  AEOFM 

TMETR«ALFAR-W  t  NCR+GAMR 
ISFG  «  5 
GO  TO  20 

C  INTEGRATION  FROM  TOUCHDOWN  TO  NOSEDOWN 

250  TtMS-0,0 
ThRTL-O.OO 
XE-0.0 
ZF-0.0 

XEOOT  -VAPF4C0S  <  GAMR I 

ZEPOT-O.O 

GAMR-0.0 

ALFAR-THETR+W ! NCR 
INOEX-1 
I  NT-2 

Q— ROTATN#DTOR 
COEF-RCOEF 
DTIME-0.2 
I LOOP- l 
VTF-XEDOT 
255  VTF1-VTF 
CALL  GINTG 
ISFG  «  6 
GO  TO  20 

260  IF  (ALFAR-tf INCR)  265.270.255 
265  DTMWINCR-ALFARJ/Q 
ALFAR-WINCR 

VNDF-  VTF-(VTF-VTF| >/DTIME#DT 
XE-XE+  ( 0 ,5#  <  VTF+VNOF ) -VWF  >  #DT 
VTF-VNDF 

tims-tims+dt 

270  TIME2-T IMS 
DIST2-XE 
ISEG  «  7 
INDEX- I 

OS-PHOS*VTF*VTF 
MACH-VTF/SOUND 
CALL  AEOFM 
THETR-0.0 
GO  TO  20 
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c 


INTEGRATE  FROM  NOSEDOWN  to  stop 


300  ! NT«  l 

0TIME«0.2 

T|MS«0«0 

XF«0.0 

2F*0.0 

THETR-O.C 

zeoot«o.o 

0«0»0 

XFDOT-VNDF 
!  LOOP- 1 
30*  VTF 1 *VTF 

IF  (TIMS+TIME2.GE.TSP)  ISP*1 

IF  ( T I MS+T I ME2*GE«TBRK )  C0EF*BC0EF 

IF  (TIMS+TIME2.GE.TREV)  IREV«1 

IFUSP.EQ.l.AND.USP.EQ.l)  GO  TO  315 

IF  (tSP.EQ.O)  GO  TO  315 

JSP*1 

WRITE  (PRINT* 1005) 

1005  FORMAT ( 1HC*13X*#LIFT  DUMPERS  DEPLOYED#) 
LINE  «  LINE+2 

315  IF  (COEF.EQ.3COEF.AND.JBRK.E0.1 )  GO  TO  320 
IF  (COEF.NE.BCOEF)  GO  TO  320 
JRRK«1 

WRITE  (PRINT* 1006) 

1006  FORMAT ( 1  HO  *  1 3X  *#5RAKES  ON*) 

LINE  *  LINE+2 

320  IF  (IREV.EQ.l . AND . JREV.EQ. 1 )  GO  TO  325 
IF( IREV.EO.O)  GO  TO  325 
JREV» 1 
NfNG*l 

WRITE  (PRINT* 1007) 

1007  FORMAT! 1  HO, 13X,#RE VERSE  THRUST  ON#) 

LINE  «  LINE+2 

325  CALL  GINTG 
I SFG  *  8 
GO  TO  20 

330  t  f  ( XEDOT-VWF )  335,340,305 
335  OT  *  ( VWF-XEOOT ) / ( VTF 1 -XEDOT ) #DT I MF 
XF  ■  XE+ ( 0 »5# ( VWF+XFDOT ) “VWF ) *DT 
TIMS«TIMS-OT 
340  VTF  «  VWF 
I SEG  «  9 
INDEX* 1 

QS*RHOS#VTF#VTF 
MACH«VTF/SOUNO 
CALL  AEQFM 
GO  TO  20 

OUTPUT  SUMMARY  OF  LANDING 


350  GAMT«GAMTR#RTOO 
GDIST*DIST2+XE 
GTI ME»T IME2+T I  MS 
TO  I  ST  »GD IST+OtSTl 
TT'imE»GTIME+TIME1 
VAPX*VAPF#FTOK 

WRITE  (PR1  T *  1 05C )  GD!ST.GTIME*GAMT.DIST1 ,TIME1 *VAPK.TOISl  *TT1M£,V 


1050  FORMAT! 1 H-.l OX «17HGR0UND 
1«F6.2*4H  SEC.4X, 13M6LI0E 
2E  »«F8.1«3H  PT.5X* 

3  10HAIR  TIME 

42«4H  KTS*//*1 1X*16HT0TAL 
SF6*2*4H  SEC*  6X.13HSTALL 
RETURN 


DISTANCE  ■ »F8* 1 «3H  FT«4X« 13HGR0UN0  TIME  ■ 
SLOPE  » *F6»2 «4H  DEG.//* 1 IX. 14HAJR  D I ST ANC 

■«F6*2*4H  SEC*  5X* 16HAPPP0ACH  SPEED  >«F6* 
DISTANCE  ■*F8.1*3H  FT.5X* 12HT0TAL  TIME  >« 
SPEED  >«F6«2*4H  KTS ) 


C  OUTPUT  SECTION  FOR  DEFAULT  NOTES 


400  WO I TE ( PR ! NT  *1051 ) 

1051  FORMAT! 1 HO. 9X.#GAMMA  IS  GT.  OR  £0.  TO  ZERO  -  RETURN  *> 

RETURN 

405  WRITE (PRINT  *  1 0521 

1052  FORMAT! 1M0«9X.#TH£TA  IS  LT.  ZERO  -  CONFIGURATION  LANDS  NOSEWMEEL  F 
HOST  -  RETURN#) 

ALFAD  «  ALFAR*RTOD 

GAMD  «  6AMR#RT0D 

TMETD  ■  TMETR*RTOO 

WRITE (PRINT. 1053)  ALFAO.GAMD*TMETD 

1053  FORMAT (1 HO *12X«#ALPHA  ■#*F7»3«#  DEG#*/* 13X*#GAMMA  «»*F7*3«*  DEG#* 
1/*13X.#THETA  ■#  .F7.3*#  DEG*) 

OFTURN 


END 


SUBROUTINE  L00K<NX.XC.NY«YC.NZ«ZC.B.B2>  look 

C  6400  FORTRAN  SUPROUTINE  WHICH  DOES  FOUR  POINT  LAGRANGI AN  LOOK 

C  INTERPOLATIONS  IN  THREE  OJRECTIONS  FOR  TWO  DEPENDENT  VARIABLES*  LOOK 

LOOK 

C  XC*  VC*  ANO  ZC  ARE  INCREASING  LISTS  OF  THREE  INDEPENDENT  VARIABLES  LOOK 
C  B  ANO  82  ARE  DEPENDENT  VARIABLE  ARRAYS  CORRESPONDING  TO  ALL  LOOK 

C  COMBINATIONS  OF  THE  INDEPENDENT  VARIABLES  BV  VARYING  FIRST  XC*  LOOK 

C  SFCONO  YC  AND  LAST  ZC.  LOOK 

C  NO  «  I  *  WILL  RETURN  A  VALUE  Of  0  AS  ANS  AND  A  VALUE  OF  B2  AS  ANS2  LOOK 

C  FOR  SPECIFIED  VALUES  XA.VA  AND  ZA  OF  XC«YC  ANO  ZC*  7APR2 

C  NO  «  2*  WILL  RETURN  A  VALUE  OF  B2  AS  ANS2  FOR  SPECIFIED  VALUES  LOOK 

C  ANS.XA  ANO  ZA  OF  B.XC  ANO  ZC.  7APR2 

C  NO  •  3,  WILL  RETURN  A  VALUE  OF  B  AS  ANS  FOR  SPECIFIED  VALUES  LOOK 

C  XA«  YA  ANO  ZA  OF  XC.YC  ANO  ZC.  7APR2 

C  NO  ■  4*  WILL  RETURN  A  VALUE  OF  YC  AS  ANS2  FOR  SPECIFIED  VALUES  LOOK 

C  ANS « XA  AND  ZA  OF  B.XC  ANO  ZC*  7APR2 

C  NDFR(l)  •  GREATER  Than  ZERO*  CALCULATE  DB2DX  AND/OR  DBDX.  LOOK 

C  N0FRI2)  »  GREATER  THAN  2ER0.  CALCULATE  DB2DV  ANO/OR  DBDV.  LOOK 

C  N0FRJ3)  ■  GREATER  THAN  ZERO.  CALCULATE  D32DZ  AND/OR  DBDZ.  LOOK 

C  THIS  SUBROUTINE  CALLS  FOR  THE  FOLLOWING  ENTRIES.  LOOK 

C  LAGRA.INTP2  LOOK 

LOOK 

DIMENSION  XL(3) *XM(20) »XN< I C ) .YL<5) .YMI20) »YNI 10) .YLL(S) *  LOOK 

I  VMM(20) ,YNN( 10) «B( 1 ) <B2( 1 ) iXC( I } «YC( 1 ) «ZC( 1 )  LOOK 

COMMON/L I ST9/ ANS . ANS2.ND • XA • Y A . Z A . I NO I C . IFLAG.NDER ( 3 ) . DBOX . OBOY  «  L l ST9 

•  DBDZ . DB2DX .D820Y . DB2DZ  LIST9 

C  This  EQUIVALENCE  STATEMENT  IS  NOT  USEO  IN  The  MILSTOL  PROGRAM 
C  EQUIVALENCE  (L.Il.tK.M)  LOOK 

DVOX(W)  •  C3*W*(2»*C2<*W*3.«C1  )  7APR2 

ASSIGN  350  TO  IVHIT 

LOOK 

I NO I C* I  LOOK 

LINX  ■  LINY  e  LINZ  •  0  LOOK 

DC  10  I ■  t  .NX  LOO*; 

IFIXA-XCt! 1)1 l.?0»|0  7 APR » 

10  CONTINUE  LOOK” 

II  •  NX-1  LOOK 

IF  •  NX  LOOK 

LtNX  •  I  LOOK 

GO  TO  25  LOOK 

11  IP( 1-21 120.12. 13  look 

l 20  11“  1  LOOK 

IF  ■  2  LOOK 

LINX  •  I  LOOK 

GO  TO  25  LOOK 

12  I  *  3  LOOK 

1 3  fr»  I-NXU7.16.J6  LOCK 

16  1  •  NX-l  LOOK 

IT  ||  «  !-?  LOCK 

t  F  «  | ♦ I  LOOK 

GO  TO  25  LOOK 

20  |FIN0E»« I l*Nf.C>GO  TO  11  LOOK 

II  •  I  LOOK 

IF  *  I  LOOK 

25  GO  T0t26.26.2S. 26). NO  LOOK 

26  J I  •  1  L0&< 


A 

§. 

V 


Si' 

$ 
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JF  ■  NV 

LOOK 

CO  TO  45 

LOOK 

28  00  30  J» | *NV 

LOOK 

IFlYA-YCIJl >31 *40.30 

TAP92 

30  CONTINUE 

LOOK 

Jl  ■  NV-1 

LOOK 

jf  ■  nv 

LOOK 

LJNV  «  | 

LOOK 

GO  TO  45 

LOOK 

31  |F(J-2>3l5«32*3S 

LOOK 

315  J|  «  1 

LOOK 

JF  »  2 

LOOK 

LINY  »  1 

LOOK 

GO  TO  45 

LOOK 

32  J  ■  3 

LOOK 

35  1 F  ( J-NY 1 37 1 30  *  36 

LOOK 

36  J  ■  NY-1 

LOOK 

3T  J|  •  J-2 

LOOK 

JF  ■  J4I 

LOOK 

60  TO  45 

LOOK 

40  |F(NOFQ(2)«NE.O)GO  TO  31 

LOOK 

JI  «  J 

LOOK 

JF  •  J 

LOOK 

45  00  60  K  «  I.NZ 

LOOK 

!FIZA~ZC«K>  >51 *60*50 

7APR2 

50  CONTINUE 

LOOK 

or  |  •  NZ-i 

LOOK 

KF  •  NZ 

LOOK 

LIN2  *  1 

LOOK 

60  TO  65 

LOOK 

51  IF tF-2 >315*52 *55 

LOOK 

513  F|  •  I 

LOOK 

KF  «  2 

LOOK 

LINZ  «  1 

LOOK 

60  TO  65 

LOOK 

52  K  *  3 

LOOK 

53  1FOC-NZ >57*56.56 

LOOK 

•6  F  ■  NZ-1 

LOOK 

57  K|  •  K-2 

LOOK 

»fF  *  K*| 

LOOK 

|F»  JFLAGI63.63.373 

LOOK 

575  |Fti«F-{|Fl»*6*l>  >63*58.59 

LOOK 

58  FI  •  IFLA6-3 

LOOK 

KF  *  IFLA6 

LOOK 

60  TO  65 

LOOK 

89  Irt*t-1FCAG>393«63«65 

LOOK 

595  FI  •  |Ft*6 

LOOK 

KF  •  1FU4673 

LOOK 

60  TO  63 

LOOK 

60  |F»NOffOi3>,NE.6>60  TO  51 

LOCK 

FI  m  K 

LOOK 

KF  *i  K 

LOOK 

63  XLM )  m  XA 

TAPS  2 

KFI 1 >  m  V* 

VAPS2 

XNJ I >  •  Z* 

TAPS* 

»F(N0€£M3> ,EO.0)6O  TO  2CC 

LOOK 

IPATm  •  l 

LOOK 

60  TO  309 

LOCK 

200  IFfN'OEPt  11, £0,0)60  TO  210 

LOOK 

!P*Th  •  2 

LOOK 

60  TO  310 

LOOK 
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210  irtNO£Q12>«£0.0>GO  TO  220 
(PATH  •  3 
60  TO  320 

220  1P(NO£PI1).£0.0,ANO»NOER(2).EO«0)60  tO  222 

221  fF<MO0INO'2).NE4O>RETURN 

222  IPATH  *  4 
GO  TO  320 

300  17  •  I 

00  1100  MaK|«KP 
12  ■  1241 
310  IX  ■  I 

00  1200  L*tl.IF 
IX  •  1X41 
320  IV  •  1 

00  1300  55  m  J  1,JP 
IV  m  !¥♦! 

irtlPATH.lT.4lG0  TO  330 
IP<M0D(N0«2).N£.C 1GO  TO  330 
J  m  KY+\ -JJ 
GO  TO  340 
330  J  •  JJ 
340  CONTINUE 

irtlPATH.FO.n  GO  TO  IWHIT* 13001 
12  •  I 

no  1400  K»K|.XF 
17  •  1241 

irt IPATH.EO.r J  50  TO  1 WH 1 T  .  t 300 1 
IX  m  1 

00  1900  ! ■  1 1 1  |F 
IX  •  1X41 

390  IJK  •  |4NX»f J-t4NV«tK-| 1 1 

GO  TOt  36C.3V0 . 390 .3001 . I PATH 
360  VHt |Vi  ■  0t |JK» 

|F(K|O.LC.2}VMMt  |V1  •  62115X1 
GO  TO  430 

3T0  VNIJ21  •  atlJXi 

irtHo.uE.aivHNt tzi  *  nat  tjxi 
GO  TO  4C3 

360  n.tlXl  ■  9115*1 

triH0.uE.2m.Lt  txi  ■  02 1 1 5X 1 

1900  XU* 1*1  •  xctll 

IFtUlHX.FQ.ClGO  TO  366 
VKtJZl  a  AUlHt2.xL.VU.XAl 
irtwo.ue.2iVHwt i2>  ■  AuiNt2.xu.vuu.XA> 

GO  TO  409 

366  irtft.M-  ,|Fino  TO  400 
VNttZl  «  YU t |X1 
VHHtl2l  •  VUUdXl 
GO  TO  409 
400  IN0  *  0 

CAUL  U4GPAtXU.VU.lH0> 

VHI121  •  YUM  1 
irtN0.CT.2lG0  To  *09 
CALU  CACOAtXU.VUL.  t*X>> 

VWN1I21  ■  VUUtU 
*09  KNt;2t  «  2CtK| 

1*00  CONTINUE 

irtuiH2.so.nico  to  4oe 

VNItVl  •  AUINt2.XN.VN.ZA) 

|FtN0*UE.2lVHN| |V>  a  AUlHt2.XN.VNN.ZA} 

GO  TO  430 


LOOK 

UOCK 

LOOK 

uoox 

LOOK 

LOOK 

UOOX 

UOOX 

uoox 

uoox 

uoox 

LOOK 

uoox 

uoox 

uoox 

uoox 

UOOX 

uoox 

uoox 

uoox 

UOOX 


uoox 
uoox 
uoox 

uoox 
uoox 
uoox 
uoox 
uoox 
uoox 
UOOX 
uoox 
UOOX 
UOOX 

uoox 
UOOX 
uoox 
uoox 

UOOX 
TAPR2 
7AP«2 

uoox 
uoox 

U©OX 
ucox 
uoox 
UOOX 

uoox 

UOOX 

uoox 
uoox 
UOCK 

uoox  "S. 

UOOX 

UOOX  t 

TAP02 

TAPR2  '  ■! 

uoox 
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*03  |F(Kt»NE*KF)GQ  TO  *20 
VMCIY)  •  VNUZ) 

IFCND.LE.21YMMC  |V>  •  VNN<|Z) 

00  tO  *30 
420  INO  •  0 

CAUL  LAGPACXN.VN.tNO) 

VHI |Y)  •  YNCI  ) 

IFCNO.GT .2 )GO  TO  430 
CALL  LAGPACXN.YNN. INO) 

VmM||Y)  ■  YNNCI) 

*30  XM(fY)  ■  VCCJ) 

IF C  |PATH*C0«4)G0  TOC  1300. 1 03, 1300.98)  .NO 
1300  CONTINUE 

tF)L!NV.EO.O)GO  TO  438 
ANS  ■  AL I N  (  2  «  XM  .  VM  ,  V A  ) 

IF! lPATn,eO,3lOBOY  •  SLI2.XM.YM) 
IFCNO.GT. 2)G0  TO  *60 
ANS2  •  ALINI2.XM.YMM.YA) 

IF( IPATH.E0.310S20Y  •  SLC2.XM.VKM) 

00  TO  *60 

*33  IF(Jf.NE.jF)GO  TO  4S0 
ANS  •  YM(tV) 

IF (N0.LC.2IANS2  •  VMM(|V) 

00  TO  *60 

*30  CALL  INTP2C3.XM.VM.YA, ANS. Cl .C2.C3.C4) 

IF <  f PATH.EQ.3 JOPOY  •  OWOXCYA) 
|F(N0.6T.2)GO  TO  460 

CALL  INTP2I3.XM. VMM.VA. ANS2.CI .C2.C3.C4) 
IF C )PATM.£0.3)0P20Y  •  OWOX(VA) 

460  IFC |PATM.3)490,221,48C 

400  DFTUQN 

*90  VLCIX)  •  ANS 

|Ff NO.LC .2 ) YLL I  IX)  ■  ANS2 
1200  XLMX)  •  XCCI  ) 

ircLtNX.ro. OlGO  TO  498 
ANS  •  ALtNC2.XL.YL.XA) 

IFC |PATm,E0.2)030X  *  SL<£.XL«VLI 
)FiN0.GT.2)G0  TO  S22 
ANS2  •  ALtNC2.XL.YLL.XA) 
trc |PATM.C0.2)0020X  ■  SLt2.XL.VLL* 

GO  TO  322 

*93  IFCII.Ng.lFiGO  TO  310 
YNt|2)  •  VLCIX) 

|FCN0»LC.2)YNNC|2)  •  VLLUX) 

GO  TO  1100 

310  CALL  INTP2i3.XL.VL.XA, ANS, Cl. C2.C3.C4) 
|Ft|PATM.C0.2)0RPX  •  OX&XCXA) 
|FCN0.GT.2)G0  to  382 

CALL  INTP2I3.XL.VLL.XA.ANS2.CI .C2.C3.C4) 
|FC |PATH,eo.2)Q920X  «  OVOXCXA) 

322  IFC |PATh.£Q.2)GO  TO  2)0 
YNC)2)  *  ANS 
YNNCI2)  .  ANS 2 
UCO  XNCI2)  •  2CCK) 

tFCLtN2.ro. 0)60  TO  S2« 

ANS  •  ALINC2.XM.VM.2A) 

W*p2  *  SlC2.xm.yM) 

1FCN0.GT.21G0  To  200 

ANS?  *  AL I N  c  2 • XM  «  VMM . 2A ) 

0P2O*  •  SL42.XM.VMM) 

GO  TO  200 


LOOK 

LOOK 

LOOK 

LOOK 

LOOK 

LOOK 

LOOK 

LOOK 

LOOK 

LOOK 

LOOK 

LOOK 

LOOK 

LOOK 

7APR2 

LOOK 

LOOK 

7APR2 

LOOK 

LOOK 

LOOK 

LOOK 

LOOK 

LOOK 

7APP2 

7APR2 

LOOK 

7APP2 

7APP2 

LOOK 

LOOK 

LOOK 

LOOK 

LOOK 

LOCK 

7APP2 

LOOK 

LOOK 

7APP2 

LOOK 

LOOK 

LOOK 

LOOK 

LOOK 

LOOK 

7APP2 

7APP2 

LOOK 

7AP»2 

7APP2 

LOOK 

LOOK 

LOOK 

LOOK 

LOOK 

7APP2 

LOO* 

LOCK 

7APU2 

LOCK 

LOOK 


1-28 


520  IFlKt.NE.KFtGO  TO  540 
4NS  •  VNt 1Z) 

ANS2  •  VNNItZ) 

540  CALL  INTP2(5.XM,YM,ZA.ANS.C1.C2.C3.C4) 
p^z  •  owox(ZA) 
triW0.GT.2tG0  TO  200 

CALL  INTP2(5.XM,yM«.ZA.ANS2,C1.C2.C3.C4) 
PP202  •  OMOXIZA) 

60  TO  200 

06  VMH( |v j  a  XM< JVt 
100  ir J J-NY) I C | .106. »06 
tot  tr»J-l )tC7«t07.i02 
t02  1F» JJ-3tl300.l30C.104 
104  trCANS-YMCtr-l ))|303.1tS.llC 

106  !F(ANS-YM(|Vttl3CC.lt6.114 

107  irtANS«YM(tY)>113. 116*110 

110  CALL  INTP21 IV.YM.YMM.ANS.ANS2, Cl .C2.C3.C4) 
|Y  •  |V-4 

PO  til  | ■ 1 .4 
IY  •  IV41 
J  •  NV42-1Y 

111  VMM I  0 Y |  •  YC1J) 

CALL  1NTP2(|V.YM.VHM.ANS.VA.C1 .C2.C3.C4t 
PfTUBN 

113  lNO!C«2 
60  TO  116 

114  iNOtCoS 

AN?  •  YM(IY) 

116  ANS?  «  YMMJIV, 

J  a  NY42-1Y 
YA  •  YCiJt 
PFTU»N 

11*  1 Y  a  IV-I 
60  TO  116 
FNO 


LOOK 

LOOK 

LOOK 

7APR2 

7APR2 

LOOK 

7APR2 

7APR2 

LOCK 

LOOK 

LOOK 

LOOK 

LOOK 

LOOK 

LOOK 

LOOK 

LOOK 

7APP2 

7APP2 

7APP2 

7APR2 

7APR2 

7APT32 

LOOK 

LOOK 

LOOK 

LOOK 

LOOK 

LOOK 

TAPQg 

7APP2 

LOOK 

7APH2 

7APP2 

LOOK 
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SU»StOyTJN£  NKBP3tOX.Y|N*X|N» YSAV.XSAVt 

6*08  FORTRAN  SU6*OuTJN£  TO  00  A  N£*TGN  0<  LlN£AO  CCNVf SCONCE . 

ir  tox t  to.20.io 

10  Y$*V*Y|N 
XSAV«X(N 
X!*J«X!N*0X 
PX»0.0 
Cf'  Tt>  30 

20 

XtN4X|N-Y|N*Ut|N>XS4V)/<Y  1N-YS4Y1 
YKAV.X 
Y$AV*Y1N 
30  OfTyBN 
e«i 


Aft 

Aft 

Aft 

AA 

AA 

AA 

AA 

AA 

AA 

AA 

AA 

AA 

AA 

AA 

AA 


I 

| 
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SUBROUTINE  SKIP 

6400  FORTRAN  PACE  EJECT  SUBROUTINE* 

INTEGER  PACE. PRINT, CARO 
PEAL  UFT.NACH.NU 

DIMENSION  HEAD 1 1 6 1 . NAME < I C > 


AC 

AC 


AC 


AC 


COMMON  /CSST1/  L INC .PAGE »C IMl T .CARO. PR I NT . INP. I PUNCH. I DATE .HOI 60 1  AC 
COMMON  /CIST2/  OECTO.ACFAR.THRTC.GAMR.HIND.VTF.WT. HE, OT. T|MS, KF.FU  AC 
I EC . MACM . VWF ,C IFT. DRAG . TmRST  «  RCF , EN6N0 , THV . I N  AC 

COMMON  /C I ST3/  CC.CO.S.OS.THIR.ACTHR.CCKAX.AR.CCAR.CMACG.ACPHO.CX.  AC 
1CZ.0. lMOM.ttINCR.CMDMP.OSC.ACMXO.VSF.BW  AC 

COMMON  /CISTS/  SIG.S0UN0.NU.TEMR.PAM8 . I ATM.TEMF .dsodw.drho  AC 

data  INAMEC I) * I»l • |u»/lCHSTAN0ARD  A » 9HTM0SPMERE « I CHTROP I C  TEM.9HPE  AC 
|PATURES.IOm  POLAR  TE«.9mPERATURES.IOM  HOT  TeM.9HPERATURES.lOM  C  AC 
?Ocr*  TFM.9HPERATURES/  AC 

pfiTA  I HO (I). 1*1. 59  J / 1  OH  » 10H  AfRSPEEO.lOH  ACPHA.10H  AC 

IALT  ERROR  *  J  PH  ALTITUDE. IOH  AX8.J0H  A20.10H  CC*  AC 

CO. IOM  DISTANCE.1CM  ORAG.IOM  ELEVATOR. 10H  FUEL  U  AC 

35Pn»!0H  FCT.  PATH.JOH  FUEL.JOh  FUEL  FCOW.IOH  HEIGHT. IOM  H  AC 

400*  HT,,ICmhORIZONTAC»ICH  CIFT.lOHCOAD  FACT. .IOM  MACH  N0.*)0  AC 

SHnOR.  ACCEL* I  oh  0S.10H  RANGE.jOHR/C  ACTUAL. $0M  SEGHEN  AC 

6T.I0M  SEP. ICh  SEP/WF.Ium  SP.  RANGE ♦ IOHTAN.  ACCEC*)OH  TmR  AC 

TOTTcE.lCM  ThETA.IOh  THETA  OOT.IOHTM  DSC  OPT.IOh  TfME.IOH  AC 

e  WEIGhT.ICm  (DEGREES) . ICM  IF£ET)*ICm  (F/S/SI.IOM  (FT/CB).  AC 

91 OH  CFT/MJNJ.IOH  (FT/SIC) *  I  CM  IGi.lCH  (KNOTS). iOH  IC0/  AC 

Ahbj.jCk  (MINUTES) . ICM  <n  M J /c8 ) • t Oh  (N  MfCESI.tOM  (PERCENT )* |0h  I  AC 
»PQUNOS).|CH  THRUST, ICHRXC  INSTAN.IOM  TEMP. • l OH  (OEG  F)*lO  AC 

ah  FLAP. I  OH  ( SECONDS ) . t OhYmST  AcPhA.IOh  (OEG/SEC)/  AC 

AC 

the  input  consists  of  two  alphanumeric  title  cards 

CARD  1  -  COLUMNS  2  through  2C  ARE  RESERVED  FOR  ENGINE 
IDENTIFICATION 

COLUMNS  2t  THROUGH  ac  ARE  pointed  as  a  TITLE  CINE 
CAO0  2  •  COLUMNS  I  THROUGH  80  ARE  PRINTED  AS  A  SECOND 
TITLE  UNE 


IF  IIN.NE.II  GO  TO  1C 
OFA0  (CARD. I  50)  (KfAD* !)«!•» .16) 
130  FORMAT  i IX.A9.7AtO/eAlO) 

OfYoRN 


AL 

AC 

AC 


io  calc  sec©ko»t»  ac 

WOITE  (PRINT.) AO)  IQATt.T.PAGE  AC 

140  FORMAT  (IMl ,4X,?fNAl«Cfi*PT  PE«F0«H4NCF.*9X,JI0.I0X,FJ0.2.SX.*JMP*GE  AC 
I  .13)  **• 

hFNG*En6N0  AC 

WRITE  (PRINt.iSD  NENG.HgAOtl »*k£AD(2S .$  PRE 

ISO  FORMAT  («H  .2SX.S2*  t' n& i NS • A9 . A l 0*  S  *  4F5.04  SO  Ft  •»  PPL 

IF  (  I  ATM.EOv.f*)  GO  TO  Sf,  .  '  .AC 

IF  (5*Tm,EO.O«AiO.DT,n€,v,)  GC  TO  30  *L 

N*2»(IATM*)>  AC 

WRITE  (PRINT,  }60)  NAH£fN«i  ).NA«t*{N|  AC 

160  FORMAT  1 1HW.65X.AI0.A9)  AC 
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m 

m 


60  TO  40 

20  WPITE  {POINT.  170  TfMF 

170  FOPMAT  I lM*.e»Xl2HTEMPEPATU«E**F5*l *SH0£6  F> 
00  TO  40 

30  WPJTE  {POINT.  180  DT 

160  F09MA?  Ci«4*84K5HTEMP».FS.1.12H0E6  POOH  STCJ 
40  XPJTE  {POINT. 1901  (MEAD{ I ) «1«3.16) 

190  F OP*AT  ( 35X  «6A  J  0/25X • 8A 1 C ) 

PAGE "PAGE* 1 

IF  (JNP.6T.9J  60  TO  SC 

60  TO  (6C.60.60.60.00.90.  60.60.60).  JNP 
SO  JP-lNP-tS 
60  TO  120 
60  L|NC«4 
DfTUON 

BO  WOITE  (POJNT.21CJ 
210  F0OMAT  (20X.«ST0L  TAKEOFF#) 

60  TO  120 

90  WPITE  (POINT# 220) 

220  rOPPAT  (2CX.#ST0L  LANDING*) 

120  L|Nf*5 

ofTvan 


AL 

AL 

AL 

AL 

AL 

AL 

AL 

AL 

AL 

AL 

AL 


'4 


s 

3 


AL 

AL 

AL 

AL 

AL 


AL 

AL 

AL 


FNO 


"  :  •  •  :  •  ■$ 
■■■■.■■  .  S 

'  •  “L 

‘  ‘  .  ‘it 

'i 

•  :  ••  .  | 
.■  .  g 

V  .  '  .  .  T 

■  4 

;  FUNCTION  SLM.X.V)  J 

*UnT»}).-.  .-••••-  •  •,  1 

..  •  • 

5L  •  «v(  {♦)  j-y;|  }?  >  { 

orrvsj<6  .■  3 

-  •  -  •  '  'I 

Fn0 
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SUBROUTINE  SMLT2(NSCL*CN| 

C  GFMt.RAL  SIMULTANEOUS  EQUATION  SOLUTION  SOAR  OUT  INC. 

NSOi.  •  MO*  OF  SOLUTIONS!***  61*  CN  •  COEFFICIENTS*  EQUATION  FORM 
CNIItll  •  CN( I *2lAM IFCNI I *3)A(2)a..«CN( |*NSOL+t lA(NSOL) 

WHERE  I  TAXES  VALUES  t  THROUGH  NSOL.  ANO  A<|>  ARE  SOLUTIONS. 
SOLUTIONS  arc  RETURNED  AS  CNIItU* 

DIMENSION  CNI6.T) 


1 1 -NSOL 
CO  TO  20 
10  ll-IL 
20  JJ*]|*| 

DO  60  1 •  1  •  I  f 
IF  CNM.JJI.EO.0.0 

IF  I.GT.I  A DO  £0  I  TO  EO  |* 

IF  I *E0«  t  FIND  LOWEST  ||*N)  WHERE  CNt  1 * JJI.NE.O.O. 
AND  AOO  £0  N  TO  EO  S. 

IFICNC I *JJ> .NE.O.OlGO  TO  30 
DO  32  N* | * | | 

|FICN<N.JJ|.N£.0*0|G0  TO  36 
32  CONTINUE 

WRITE (6* I COOIJJ 

JOOO  FORMAT I I0X6mC0LUMN«2X* I2.2X6NALL  ZERO! 

STOP 

36  oo  ye  m«i %jj 

36  CN? I *M|  ■  CNI I <H)*CNIN*M) 

DIVIDE  ROW  ©V  CNM.UUI.  LEAVING  I*  IN  LAST  OR  JJ  COLUMN  OF 
CURRENT  matrix. 

TO  OCPal  ./CNM  *JU> 

DO  AO  J*| i JJ 
40  CM(|*J.*CN(I.U|«RCP 
60  CONTINUE 

STORE  ECUATIQN  FO «  CALCULATING  SOLUTION  JN  EMPTY  COLUMN. 
LAST  NUMBER  STORED  IS  F|HST  SOLUTION. 

DO  63  |*t. (I 
63  CN« | » JJi  «  Qni i , | | 

IL*1 1-1 

IF  ML  I  60*00. TO 


L 

L 

L 

L 

L 


L 


L 

L 

L 

L 

L 

L 

L 

L 


L 

L 


¥MFN  IL  A  0.  YM*  FIRST  SOLUTION  IS  COMPLETE. 


ASSAM' l**0  t  IN  EAC«  LAST  COLUMN,  SUBTRACTIN',  SUCCESSIVE  EQUATIONS 
RESULTS  IN  ONE  LESS  EQUATION  EACm  W|T-  Cn£  LESS  UNKNOWN. 

TO  It  BO  I • I « IL 
12*5*1 

DO  f}0  J*l.|  | 

60  CnM  »  Ji*C««  f  •  |T.  Ji 

GO  tO-iO 

60  IF  inSCl-M  500*120*100 
C  CALCULATE  remaining  SOLUTIONS*  store  IN  FIRS?  COLUMN. 

*06  CO  ISO  J*2.NSOL 
J|T|*I 

CNt  S  < I l»CNt !.J|| 

ne  s to  u*i«jf 

; 

? |C  FNf|, 8  *»CN*  a .) |«CN,U2«U| |*CNC u*5  j 

*20  wtijtoi 
TNO 


L 

L 

L 

L 

L 

L 

L 

L 

L 

L 

L 

L 

L 

L 

L 

L 

L 
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SUBROUTINE  TAKEOFF 

6400  FORTRAN  IV  TAKEOFF  SUBROUTINE 

STOL  TACTICAL  AIRCRAFT  INVESTIGATION  GROUND  RULES  -  JULY  1972 

THF  BOUNDARIES  OF  THIS  PROGRAM  ARE  - 
VR  GT  OR  EQ  TO  VI  AND  VMC  GROUND 
VLO  GT  OR  EQ  TO  VSC  *  VMC  AIR 

ALPHA  AT  LIFTOFF  LT  OR  EO  TO  GALMX  (AERO  INPUT) 

ALPHA  DURING  ROTATION  LT  OR  EQ  TO  ALPMX  (AERO  INPUT) 

THIS  PROGRAM  CALLS  THE  FOLLOWING  SUBROUTINES* 

ATMOS*  SK IP*  HEAD*  THST 1 «NWRP2  « AERO  1 • AEQFM*G INTG 

INTEGER  CARD* PR  I NT* PAGE 
REAL  L!FT«KTOF,MACH«NU*NMtOF 

DIMENSION  NL(B) 

COMMON  /LIST1/  L INE  «PAGE  *L IM I T  * CARD  * PR  I  NT  t I NP« I PUNCH* I  DATE  «  HD (60 )  F 
COMMON  /L I ST2/  OELTD  * ALFAR . THRTL . GAMR . H I ND ♦ VTF  t  WT  t  HF ♦ DT , T I MS «  XF ♦ FU 
1 FL « MACH . VWF  «L I FT ♦ DRAG ♦ THRST  *  RCF *  ENGNO • THV* IN 
COMMON  /L I ST3/  CL . CD . S . QS ♦ TH I R ♦ ALTHR • CLMAX « AR . CLAR * CMACG . ALPHD *CX * 

1 CZ.Q « I MOM  »  W I NCR , CMDMP  *  QSC . ALMXR ♦ VSF . SW 
COMMON  /LIST4/  FFS* JPOW  *TGROS  *DWA  » THMOM  *  THREQ 
COMMON  /LISTS/  S I G  *  SOUND  *  NU  *TEMR«PAMB« 1  ATM »  TEMF  *  DSODH  «  DRHO 
COMMON  /LIST6/  RTOD.DTOR*KTOF.FTOK*NMTOF.FTONM,RHOZ*RHOZ2*GZ«PZ.TZ 
COMMON  /LISTS/  I NDEX  *COEF  »GN*GT*AN*AT«AX«AZ* I TRM  * VUPPER 
COMMON  /L I  ST 1 5/  I ALMX , 1 NT , XEDOT *  ZEDOT . XE  *  ZE , I  LOOP . DT 1  ME , I  AT , FL I FT . 
1DFLP.C4HT.GALMXR 

COMMON  /lgeom/  shsw*scr«Qdot  *RTH.DTHTH*RWA*DTHWA  *DSTHZ  *dwlhz* i yy  *x 

icg*zcg*cbar.pfn, si nth .costh. angle *sinal»cosal*rhz»dthhz*vtsq*udqt. 

2WD0T  *THETR*U«W 

COMMON/CONTROL/ J.F  I G  •  I  REV  *  I  SP « NENG 
COMMON  /LIST99/  I ERR 

DATA  VSC*OFLP*RCOEF,8COEF,ROTATN*ROTPT*TIMR,TIMB*DGLO/1 .05*25. *. 1 ♦ 
IC«3*8.*0»* l.*2.*0.1C/ 

DATA  (NL( I ) *1«I .e)/l0H0556370226*10H3323082052.IOH3938514542,l0H33 
1 4 DO 15151  * 10H3610171403, 1CH34  31091 158* 10H5739393838. 10H594001513B/ 

OFF  1 N I T I  ON  OF  VARIABLES  IN  NAMELIST  • TAKEOr 1  • 

VCS  ■  RATIO  OF  LIFTOFF  SPEED  TO  AIR  MINIMUM  CONTROL  SPEED 
DGLO« INCREMENTAL  LOAD  FACTOR  REQUIRED  AT  LIFTOFF 
OFLP-FLAP  SETTING  IN  DEGREES 
RCOEF  »  ROLLING  COEFFICIENT  OF  FRICTION. 

BCOEF  «  COEFFICIENT  OF  BRAKING  FRICTION 

VMCGKa  MINIMUM  CONTROL  SPEED  ON  THE  GROUND  (KNOTS) 

VMCAK  •  MINIMUM  CONTROL  SPEED  IN  THE  AIR  (KNOTS) 

RCTATN  ROTATION  RATE  (DEGREES  PER  SECOND) 

T  t  mf  «  REACTION  T I  MR  FOR  ENGJNF  FAILURE  (SEC) 

T I  MB  «  BRAKING  DELAY  AFTER  «TIMR«  (SEC) 

ROTPT  EQUAL  TO  ZERO  SUPRESSES  PRINTING  OUTPUT  FOR 
SEGMENTS  4*6  AND  8. 

THF  FOLLOWING  VALUES  ARE  ENTERED  AT  TIME  OF  LOADING  AND  ARE  USED 
UNTIL  OVERRIDEN  BY  READING  THE  APPROPRIATE  VARIABLES  IN  TAKEOF’. 
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VSC-1.05 

dglo=o«  jo 

OFLP-25. 

PC0EF*0# 1 0 

0COEF.O.3O 

TIMR*1,0 

TtMB«2,0 

ROTATN*8,0 

ROTPT  *  l  • 

NAMEL I ST/TAKEOF 1 /VSC . DFLP . RC  OEF . BCOEF *  VMCGK  «  VMOAK  .ROTATN « 
jOOTPT*  f  I  MR  »  T  t  MB i DGLO 

IF  ( IN.NE.l )  GO  TO  10 
RFAD  ( CARD . TAKEOF 1 ) 

WRITE  (PRINT, TAKEOF1  ) 
ofturn 

10  HFRUN»HF 

MF«HFRUN+C4HT 
IBAL-0 
CALL  ATMOS 
PHOS*RHOZ2*S 1 G*S 
VMCA  «  VMCAK*KTOF 
VMCG  *  VMCGK*KTOF 
V1F«VMCG 
VI* VMCGK 

WRITE  (PRINT, 914) 

9 1  A  FORMAT! 1  HI «3X.*0UTPUT  DEFINITIONS  -  TAKEOFF*. /,4X,*S£GMENT 
1*1  DETERMINATION  of  liftoff  angle  of  ATTACK*,/, IAX, 
INCREASE  VLO  TO  AVOID  GROUND  CONTACT*./, IAX, 

R/C  AVAILABLE  AT  LIFTOFF  -  1  ENGINE  OUT*,/, IAX, 

VLO  TO  VR  INTEGRATION  STEPS* ,/, 1  AX, 

CONDITIONS  AFTER  VLO  TO  VR  I NTEGRAT ION*,/, 1  AX , 

VR  TO  VF  INTEGRATION  STEPS*, /, 1 4X , 

CONDITIONS  AFTER  VR  TO  VF  INTEGRAT I  ON* ,/. I AX. 

VF  TO  STOP  INTEGRATION  STEPS*,/. I  AX. 

CONDITIONS  AFTER  VF  TO  STOP  INTEGRATI ON* ./, 1  AX. 

V*0  TO  VR  INTEGRATION  STEPS* ♦/, I  AX, 

CONDITIONS  AFTER  GROUND  RUN  TO  VR* ) 

LINE-LIMIT+1 
GO  TO  ISO 


PRE 


*. 


2*2 

3*3 

4*A 

5*5 
6*6 
7*7 
8*0 
9*9 
1*10 
2*1  1 


OUTPUT  BLOCK, 

SO  IF  (ROTPT. NE. 0, )  GO  TO  30 

IF  ( ( ISEG.EO.A) *OR. ( ISEG.E0.6I.0R, ( I SEG.E0.8 ) .OR. ( ISEG.EQ. 10) > 
1  GO  TO  60 

30  IF  (LINE-LIMIT)  30, AO. AO 
AO  CALL  SKIP 

CALL  HEAD (NL ) 

WRITE  (PRINT, 913)  HD(27( 

913  FORMAT  ( 1 H+ , 1 X , A  1 0 ) 

50  VTK«VTF*FTOK 
GAMD«GAMR*RTOD 
RCM*60,*VTF*SIN(GAMR) 

THETD«THETR*RTOO 

THDTD*Q*RTO't 

4LFAD*ALFAR*RT0D 

althd«althr»rtod 
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WRITE  (PRINT. 91 2 >  I SEG. HF.OFLP. WT . VTK .RCM.THETO. AN.CL .LIFT . TGROS.T 
1 IMS.XE.ZE. GAMO .  ALF AO . THOTO .  AT . CD . DRAG .  ALTHD 
L|NE»L|NE+3 

912  FORMAT  ( 1  HO. I  1 1 ,6F 1 2.3. 2F 12,5, 2F 12.3/12X.5F12, 3.3F12.5 .2F12 .3) 

60  GO  TO  (250, 210. 570, 610. 67C.68C, 700, 730.750. 820. 950), ISEG 

C  SET  REQUIRED  LIFT-OFF  SPEED  (POWER  ON  STALL  LESS  1  ENGINE) 

160  |N»0 

CALL  AERO  1 

ALTMR  «  ALPHD-W I NCR+TH I R 
VTF |  ■  VSF 
ENGNO  «  ENGNO-1 . 

NFNG«2 
THRTL-1 .00 
MACH  »  VTF 1 /SOUND 
VTF  ■  VTF1 
NSTEP  «  0 

170  CALL  THST1 
CALL  AERO I 

IF ( JF I G.EQ.2 . AND . ( VTF+VWF ) ,GT« 1 . ) CLMAXpCLMAX- ( ( T GROSES  IN ( ALTMR ) )/ ( 
1VTF*VTF*RH0S) ) 

VTF  »SQRT ( (WT-TGROS*SIN(ALTMR> )/<RHOS#CLMAX) ) 

YVTF  P  (VTF I —VTF ) /VTF 

!F( ABS ( YVTF ) -0 .000 1 )  174.174.171 

171  IF(NSTEP)  173.172.173 

172  DMAC  ■  MACH— VTF/ SOUND 

173  CALL  NWRP2(DMAC. YVTF, MACH, YSAV.XSAV) 

NSTEP  »  NSTEP+1 
IF(NSTEP.GT.iO)  GO  TO  161 

VTF 1  *  MACH4S0UN0 
VTF  «  VTF1 
GO  TO  170 

161  VTF  a  MACH4SOUND 

174  VSF  ■  VTF 
IN«2 

VLOF  «  VSC^VMCA 

IVSC«0 

!LOF«0 

VTF»VLOF 

VSK«V5F*FT0K 

IF (VTF »GE* VSF* SORT ( 1 ,+DGLO) )  GO  TO  25 
ILOF« 1 

VTF«VSF*SQRT  < 1 ,+DGLO) 

25  VLOF«VTF 

VLOK»VLOF»FTOK 

C  DETERMINE  ANGLE  OF  ATTACK  AT  LIFTOFF, 

190  VTF.VLOF 
THRTL«1 .00 
ISP  «  0 
I  REV  «  0 
COEF  »  RCOEF 
TIMS*0.0 
XF«0.0 
Zf-0.0 
0  «  0,0 
IFRR.O 
I  AT«0 
t  ALMXpO 
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INTEG*0 

HF-HFRUN+C4HT 

GAMR»0, 

alfar*wincr 

IN0EX«5 

OS»RHOS*VTF*VTF 

mach-vtf/souno 

GN«1.0 

GT«0.0 

tall  afofm 

200  ALFLO»ALFAR 
THETR* ALFAR-W ! NCR 

!F  ( ALFLO— W I  NCR  1  900.202*201 

201  IF(ALFAR.LE.GALMXR)  I SEG* 1 
IF(ALFAR.GT.GALMXR)  ISEG«2 
GO  TO  20 

C  ADJUST  VLOF  TO  AVOID  GROUND  CONTACT 

210  IVSC  »  IVSC+1 

1F( ILOF.GT.O)  GO  TO  21 1 
VSC1  »  VSC+FLOAT<  IVSCl/lOO. 

VTF«VMCA*VSC1 

WRITE  (PRINT. 910)  ALFAO 

WP I TE ( PR  I NT .917)  VSC1 

917  FORMAT ( 1H  .9X.16HVL0  INCREASED  TO .F7.3 *2X . 1 2HT I MES  VMCAIR) 
LINE  ■  LINE+3 
GO  TO  25 

211  VTF.VSF# < SORT ( l ,+DGLO) +FLOAT < IVSC 1/1 00.  ) 

VSC1 “VTF/VSF 

WR I TE ( PR I  NT .910)  ALFAD 
WR I TE  <  PR  I  NT .911)  VSC1 

911  FORMAT t 1 H  .9X.16HVL0  INCREASED  T0.F7.3.2X. 12HTIMES  VSTALL) 
LINF  ■  LINE+3 
GO  TO  25 


C  DFTFRMINE  R/C  AVAILABLE  at  LIFTOFF, 

250  t  N0EX«3 

ALFAR^WINCR 
0S«RH0S*VTF*VTF 
MACH«VTF/S0UND 
CALL  AEQFM 

THETR« ALFAR-W I NCR+GAMR 
I SEG*3 
60  TO  20 

C  INTEGRATION  FROM  VLO  TO  VR, 


570  IF  (INTEG)  590,590,580 
580  SEG2-0.0 
T | ME2*0 .0 
VRFaVLOF 
GO  TO  660 
590  HFchFRUN+C4HT 
CALL  ATMOS 
RHOS»PHOZ2*SrG#S 
T JM$«0,0 
XP«0.0 
ZF-0.0 
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I ALMX«0 
IAT.0 

ThETR»ALFLO-WINCR 

xedot«vlof 

zfdot-o.o 

INineXM 

INT«2 

0«P0TATN#DT0R 
DTIME—0.2 
I  LOOP*  1 
600  VTF 1 -VTF 
CALL.  GJNTG 
ISEG»4 
GO  TO  20 

C  CONDITIONS  AFTER  VLO  TO  VR  INTEGRATION 

610  IF  <  ALFAR-W I  NCR )  620,650.600 
620  0T«-( ALFAR-WINCR)/Q 
ALFAR«  WtNCR 
THFTR  ■  0.00 

VRF-VTF+OT#  <  VTF-VTF 1 ) /OT I  ME 
IF( VRF-V1F)  625.640,640 
625  OVLOF  «  (V1F  -  VRF) 

IF ( DVLOF-. 1 )  640.640.630 


C  EMPIRICAL  FACTOR  ADDED  TO  IMPROVE  ESTIMATE  OF  DV  NEEDED 

630  VLOF-VLOF+  .75*1 VtF  -VRF) 

GO  TO  190 

640  XF=XE+ ( .5#  <  VTF+VRF ) -VWF ) *DT 
VTF.VRF 
TIMS-TJMS+OT 
650  SFG2--XE 

T I ME2*-T IMS 
660  VTK2«VL0F#FT0K 

IF ( ABS ( VTF- VI F ) »LE .1.0)  ISEG=7 
IF ( ABS ( VTF- VI F ) »GT .1*0)  ISEG«5 
INDEX=1 

GS«RHOS*VTF*VTF 
MACH»VTF/SOUND 
CALL  AEQFM 
GO  TO  20 


C  INTEGRATION  FROM  VR  TO  VF 

670  SFG3«0,00 
T IME3*0.00 
T IMS*0.00 
I  NT  •  1 

OT I ME»-0.2 
XF«0.00 
ZF«0.00 
ThFTR»0,0 
XFPOT«VRF 
ZFOOT«O«O0 
tLOOP«I 
VTF»VRF 
675  VTF l ■ VTF 
CALL  GINTG 


PRE 


mi 

m 


I 

1 

i 

i 


•M 


ft 

sfj 

■5 


-.•■it 

$ 

'5 

’•'i 

\ 


1-37 


ISFG»6 
GO  TO  20 

C  CONDITIONS  AFTER  VR  TO  VF  INTEGRATION 

680  IF ( VTF-V1 F )  685.690*675 
685  OT  ■  DTIME#< < VTF-V1 F) /< VTF-VTF 1  I  I 
TIMS  *  TIMS-OT 

XF  »  XE-( ( VTF+V1F ) #0 «5-VWF ) #OT 
VTF  *V1  F 
690  SEG3«-XE 

TIME3— TIMS 
V1*V1F#FT0K 
ISFG*7 
INOFX«l 

OS*RHOS#VTF#VTF 
M4CN«VTF/S0UN0 
CALL  AFOFM 
GO  TO  20 

c  integration  from  vf  to  stop 

700  TIMS»C,0 
I  NT  •  1 
SFG5«0,0 
DTIME*0.5 
TIME5«0.00 
XF*0.0 
7F*0.0 
THFTR«0.0 
?Fn0T»0.0 
0«0.0 
XE00T«V1F 
I LOOP* l 
VTF-V1F 
710  VTF 1  * VTF 

IFITIMS+O.l *0T I  ME  »GE.  T I  MR )  THRTL*0*0 

I F  (  T  I  MS+C  «  1  #DT  I  ME  «LT»T I  MR+T I  MB  !  GO  TO  720 

C0FF*BC0EF 

I RFV* 1 

NFNG  *  2 

!*P«l 

720  CALL  G I NTG 
ISPG*3 
GO  TO  20 

C  CONDITIONS  AFTER  VF  TO  STOP  INTEGRATION 

730  IF(XE00T-VWF>735. 740.710 
735  DT  «  ( VWF-XEDOT ) / ( VTF1 -XEDOT ) #DT I  ME 
XF  ■  XE+ ( C .5* ( VWF+XEOOT  J-VWF ) #DT 
TJMS  ■  TIMS-OT 
740  VTF*VWF 
1 RFG*9 
INPEX-1 

QS*RHOS#VTF#VTF 
MACH«VTF/SOUND 
SFG5-XF 
TIM£5«T!MS 
CALL  AEOFM 
GO  TO  20 


■i 
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C  MODIFY  VI  TO  BALANCE  ISEG-3  ♦  1SEG-2  WITH  ISEG-5 

750  IF  <SEG2+SEG3.Le.SEG5.AND.tBAL.EQ.0>  GO  TO  795 
IF<ABS<S£G2#SEG3-SEG5).LE«  2.)  GO  TO  795 
IF( IBAL.GT.20)  GO  TO  795 
IF< IBAL.GT.O)  GO  TO  755 
V0LD-V1F 
D1-SEG2+SEG3 
P2-SEG5 
VNFW-V1F+10. 

03«0 i # <  VLOF-VNEW ) / ( VLOF-VOLD  > 

P4-02*VNEW/V0LD 
GO  TO  760 
755  D3-SEG2+SEG3 
D4«SEG5 
VNEW-V1F 
760  IBAL-IBAL41 

VtF-V0L0+<  <02-01  lAtVNEW-VOLO/IDS+DE-Dl-DA) ) 

VJ  «  V1F*FT0K 

IF< IBAL.EO.J >  GO  TO  765 

VOLD-VNEW 

Dl-03 

02-04 

765  WRITE (PR I  NT .915)  VI 

915  FORMAT ( 1 h0» 1  OX  *#F I  ELD  NOT  BALANCED# .5X . #NEW  VI  SPEED  ■  #,F7.2. 
1#  KTS#) 

GO  TO  190 

795  WR 1 TE ( PR I  NT .916)  IBAL 

916  FORMAT (1H0.10X. #F 1ELD  BALANCED  IN  #.12.#  ITERATIONS#) 

C  INTEGRATION  FROM  V«0  TO  VF 

800  DTIME-l.O 
INOEX-I 
HF-HFRUN+C4HT 
CALL  ATMOS 
PH0?»RH0Z2#SIG#S 

tims-o.o 

XF-0.0 

ZF-0.0 

tmftr-o.0 
YhRTL- l . oc 
ENGNO  ■  ENGN0+1. 

NFNG-0 

COEF  •  RCOEF 
JRP-0 
I RTV-0 

xpdot-vwf 
ZFOOT-O.O 
0-0.0 
I LOOP- 1 
810  VTF1-VTF 
CALL  GINTG 
I SrG- 1 0 
GO  TO  20 

C  CONDITIONS  AFTER  V-0  TO  VF  INTEGRATION 

820  IFIVTF-V1F)  810,840.830 

830  OT--OT JME# ( VTF-V1F ) /( VTF-VTF1  ) 


4 
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xr«  XE+ ( <  VTF+V IF)#. 5-VWF ) *0T 
TJMS*TlMS*OT 
VTF*V IF 
640  SFG1-XE 

TIMEI *T IMS 
VTK1«V|F*FTCK 
!SF6«lt 
inoexm 

OS«RHOS*VTF*VTF 
MACM-VTF/SOUND 
CALL  AEOFM 
GO  TO  20 

C  OUTPUT  SECTION  FOR  DEFAULT  NOTES 

900  ALFLO  ■  ALFLO  *  RTOD 

WRITE  {PRINT. 910)  ALFLO 

910  FORMAT ( I  HO  »  9X  ,  *A  I RPLANE  LIFTS  OFF  AT  ALPHA  »«F7.3.*  DEG*> 

RETURN 

C  OUTPUT  SUMMARY  OF  TAKEOFF 

950  GntST»SFGl*SFG2*SEG3 

GT1ME«TIME1+TIMF2+T|ME3 
AP I ST  »SEG1 +SEG5 
ATIME«T IME1+TJME5 
VLOK«VLOF*FTOK 
VRK«VRF*FT0K 
V1K*VIF*FT0K 

WRITE (PR I  NT .955 1  V1K  »VRK  « VLOK  *VSK 

955  FORMAT ( 1  HO « 1  OX «*ENG1 NE  FAILURE  SPEED  ■».F7.2«*  KTS*./. 11X.*R0TAT|0 
IN  SPEED  »*.F7.2,*  KTS*./.  UX.*LtFTOFF  SPEEO  «*.F7.2.*  KTS*«/.11X,* 
2STALL  SPEED  «».F7.2,*  KTS  (WITH  ONE  ENGINE  OUT)*) 

WPITE(PRTNT,956)  6DIST.GTIME 

956  FORMAT!  1H0.!0X,*V»0  TO  LIFTOFF*,/.  UX,*DISTANCE  **.F7.2.*  FT*,10X. 
1 *T ) ME  ■* «F7*2 . *  SEC*) 

WP | TE ( PR l NT . 957 )  AD  t  ST , AT |MF 

957  FORMAT U HO, 1  OX, #AB0RT£0  TAKEOFF* ./, II X » *D l STANCE  ■«.F7.2,*  FT*«IOX 
l,#T|ME  «*»F7»2  »*  SEC*) 

RETURN 

END 
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SUBROUTINE  THSTl 

6400  FORTRAN  SUBROUTINE  TO  CALCULATE  THRUST  AS  A  FUNCTION  OF 
TRUE  AIRSPEED  IN  KNOTS 

NOTE,.,  THIS  PROGRAM  USES  THE  COMMON  LIST  ELEMENT  -THMOM-  (LIST4)  AS  THE 
GROSS  THRUST  FROM  THE  BLOWING  SLOT 

Integer  card* print, page 

REAL  L!FT»MACH,KTOF,NMTOF«!YV 

common  /L IST1 /L I NE.PAGE  *L !M I T  *CARD  » PR I  NT . INP* IPUNCH. I DATE  * HD ( 60 > 

COMMON  /L I ST2/DELTD  «  ALF AR  *  THRTL . GAMR » H 1 ND  *  VTF  «  WT , HF , DT . T 1 MS  *  XF  » 

*  FUEL*MACH*VWF.LIFT*ORAG»THRST,RCFoENGNO,THV* in 

COMMON  /L|ST3/CL*CD*S»0S«THIR*ALThR*CLMAX*AR»CLAR*CMACG*ALPHD*CX* 

*  CZ*0* IM0M.WINCR*CMDMP,0SC«ALMXR,VSF»BW 
COMMON  /L I ST4/FFS *  UPOW  «  TGROS  *  DW A  *  THMOM *  ThREO 

COMMON  /lists/  sig*sound*nu,temr,pamb, iatm,temf.dsodh,drho 

COMMON  /L I ST6/RT0D . 0T0R , K TOF . F T0K , NMTOF , FTONM , RHOZ  » RH0Z2 , GZ  *  PZ  *  TZ 
COMMON/LGEOM/SHSW  *  SCR »  QOOT ♦ RTH, DTHTH ♦ RW A . DThWA  * DSTHZ  tDWLHZ . I YY * 

1  XCG*  ZCG*CBAR  *PFN,S INTH*COSTH  * ANGLE  *S INAL  *COSAL  *RHZ*DTHHZ» VTSQ* 

2  unoT ,WOOT »THETR*U,W 

COMMON/CONTROL/JFIG* i rev* isp*neng 

DIMENSION  VKS<20) *THST(20> ,RORG(20) *T!0L(20) »TSLT(20>  *TREV(20) • 

I  DWMEI20) 

DATA  ENGNO. SCALE. NENG/4.U, 1.0*1/ 

DEFINITION  OF  VARIABLES  IN  NAMELIST  »THTl* 

N-  NUMBER  OF  X.Y  POINTS  IN  EACH  TABLE 

VMS-  THE  VELOCITY  TABLE  FOR  THE  PROPULSION  TABLES  IN  KTAS 
(USED  AS  THE  INDEPENDENT  VARIABLE  IN  ALL  TABLES) 

TMCT-GROSS  THRUST  TABLE  AT  MAX  POWER  (IN  L8S»> 

RDRG-RAM  DRAG  TABLE  *T  MAX  POWER  UN  LBS.) 

TIDL-6R0SS  THRUST  TABLE  AT  IOLE  POWER  (IN  LBS.) 

TSLT-GROSS  THRUST  AT  THE  SLOT  EXIT  (IBF  CONFIGS)  AT  MAX  POWER  (LBS) 
TPFV-MAX  REVERSE  THRUST  (EXPRESSEO  AS  A  NEGATIVE  VALUE  -  IN  LBS.) 
OWME-WINDMILLING  DRAG  for  a  OEAO  ENGINE  (IN  LIS. ) 

TH I O-THRUST  VECTOR  INCIDENCE  REF.  TO  A  WATER  LINE  IN  DEG. 

ENGNO- THE  NUMBER  OF  ENGINES 

SCALE-SCALING  FACTOR  FOR  THE  PROPULSION  DATA 
NFNG»0  NO  REVERSE  THRUST 

= 1  ALL  ENGINES  REVERSING 
•2  ENGINE  OUT  REVERSING  PROCEDURE 

ThF  FOLLOWING  VARIABLES  APF  ENTERED  AT  TIME  OF  LOADING  AND  ARE  USED 
UNTIL  OVERBIDDEN  BY  READING  THE  APPROPRIATE  VARIABLES  IN  THT1 
ENGNO  •  4.0 
SCALE  ■  1.0 
NENG  ■  l 

NflM£L 1 ST/TWT 1 /N . VK S  * ThST . RDRG « T l DL * TSL T * TR£  V . Th I D ♦ ENGNO .SCALE . NENG 
1  *DWMC 


!F( IN-1 )5, 1 *5 
l  OEAO (C ABO «ThT1 ) 
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WPITCtPPINT.THTJi 

THIR*THJD#0T0R 

ENG*ENGNO 

»FTU»N 

S  VTK*VTF*FTOK 

CALL  PIND<N«VTX*DWA»VKS«RDRG> 

CALL  FJN0<N.VTK«TMRST»VKS.THST > 

CALL  FJND<N.VTK,wMO*VXS«DWM£ ) 
IFUREV.GT.O.)  60  TO  SO 
THMOM  ■  o»oo 

fF< JFJ6.LE.2)  60  TO  10 

CALL  F  J  NO IN  « VTK  < THMQM  < VKS • TSLT ) 

10  IFfTMRTL.EQ.l «0)  GO  TO  25 

CALL  FJNOCN.VTX.THJDL.VXS.TIDL) 

RATIO*  (  TMIDL+THRTLa ITHRST-TH10L  )  1/THflST 

THRST*THRST#RATJ0 

0WA*0WA#RAT|0 

tnmom*thmom*rat JO 

2S  TGROS»THPST#ENGNO#SCALE 
OWA*OWA»ENGNO#SCALE 
1FINFNG.F0.2J  owa«dwa*wmd*scalp 
T HMOM* THMOM «CNGNO*SC ALE 


RETURN 

SO  CALL  F|ND1N.VTK«THREV»VKS,TREV) 

CALL  FINO(N.VTK,THJDL.VXS.TIOL> 
IF<HENG.EO.O)GO  TO  63 
IFJNENG.EO.l )G0  TO  60 
ThRST*THPEV*(2,0*FL0AT  < IFIXIENGN0/2 
75  TGROS*THPST*SCALE 

OWA«OW A  «SC ALE  *ENGNO 
IFJNENG.EQ.2)  DWA«OWA+WMO*SC ALE 
TmMOMbO.O 
rftupn 


.♦*0001 ) > )+THJiX 


60  THPST*TMqev»{EMGN0-2.)-*.THi0L#2. 

GO  TO  73 

65  0vA*A6S(0WA*TH JOL/THRST ) 
ThRST*Thjdl*EN6n0 
GO  TO  75 
END 
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SUBROUTINE  TRIMIN.A.V.PV.T.K)  K 

FORTRAN  SUBROUTINE  WHICH  TRIMS  N  VALUES  OF  ACCELERATION! A I  TO  ZERO  K 
OV  VARYING  N  INDEPENDENT  VARIABLES  IV).  K 

K 

PV  «  IS  A  SET  OF  PERTURBATION  INCREMENTS  FOR  IV)  TO  USE  FOR  K 

ESTABLISHING  'DERIVATIVES * •  K 

T  IS  A  SET  OF  TOLERANCES  TO  THE  VALUES  OF  I A)  WHICH  MUST  K 

BE  SATISFIED.  K 

THIS  SUBROUTINE  CALLS  SUBROUTINE  SMLY2.  K 

K 

DIMENSION  AIN)«VIN)«PV!N).TIN)t0lGt7)*DVI6)  K 


EQUIVALENCE  <D{ 1 • i ) «OV< 1 ) ) 


IF  IK)  10*10*60 
10  00  20  t  a  1 «N 

IF  <  ABS I A 1 1 ) ) — T It))  20*20*30 
20  CONTINUE 
KaNAl 
RETURN 
30  Ka| 

DO  AO  t  a | *N 
40  011*1) a»A  <|) 

SO  VIK)aV|K)«PVIK) 

RETURN 

60  00  70  |a).N 

70  011 *K4| ) a  I At  I ) +01 1*1) )/PVIK) 
V(K)aVIK)-PVIK) 

IF  IK-N)  80*90.90 
BO  KeK*l 
GO  TO  SO 

90  CALL  SMLT2IN«0) 

DO  100  t  a | *N 
100  V(|)aVt|)*OV(|) 

Kao 

RETURN 

END 
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APPENDIX  H 

PROGRAM  AND  SUBROUTINE  FLOW  CHARTS 


The  following  flow  charts  are  contained  In  this  appendix. 

Title 

Description 

Page 

MILSTOL 

Main  Program 

II-l 

MILSTOL 

Overlay  Dump  Program 

11-3 

AEQFM 

Equations  of  Motion  Subroutine 

II-4 

AEROl 

Aerodynamic  Data  Subroutine 

n-8 

ALIN 

Linear  Equation  Function 

U-10 

ATMOS 

Atmospheric  Properties  Subroutine 

11-11 

FIND 

1- Dimensional  Table  Lookup  Subroutine 

11-13 

GILL 

Integration  Subroutine 

n-14 

GINTG 

Integration  Driver  Subroutine 

II- 15 

HEAD 

Page  Heading  Subroutine 

II-l? 

INTP2 

Curve  Fitting  Subroutine 

11-18 

KABD 

Hyperbolic  Curve  Fit  Solution  Subroutine 

n-19 

LAGRA 

Lagranian  Interpolation  Subroutine 

II- 20 

LANDING 

Landing  Trajectory  Driver  Subroutine? 

11-21 

LOOK 

3- Dimensional  Table  Lookup  Subroutine 

u-u 

NWRP2 

Newton- Wr apson  Intimation  Subroutine 

n-3? 

SKIP 

Pago  Eject  Subroutine 

11-38 

SL 

Linear  Slope  Function 

11-40 

SMLT2 

Simultaneous  Equation  Solution  Subroutine 

n-41 

TAKEOFF 

Takeoff  Trajectory  Driver  Subroutine 

11-43 

THST1 

Propulsion  Data  Subroutine 

Ii-61 

TRIM 

Aircraft  Trimming  Subroutine 

11-53 
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